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Abstract 

The hypercentral Constituent Quark Model (hCQM) for the baryon 
structure is reviewed and its applications are systematically discussed. 

The model is based on a simple form of the quark potential, which con¬ 
tains a Coulomb-like interaction and a confinement, both expressed in 
terms of a collective space coordinate, the hyperradius. The model has 
only three free parameters, determined in order to describe the baryon 
spectrum. Once the parameters have been fixed, the model, in its non 
relativistic version, is used to predict various quantities of physical in¬ 
terest, namely the elastic nucleon form factors, the photocouplings and 
the helicity amplitudes for the electromagnetic excitation of the baryon 
resonances. In particular, the dependence of the helicity amplitude is 
quite well reproduced, thanks to the Coulomb-like interaction. The model 
is reformulated in a relativistic version by means of the Point Form hamil- 
ton dynamics. While the inclusion of relativity does not alter the results 
for the helicity amplitudes, a good description of the nucleon elastic form 
factors is obtained. 
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1 Introduction 

The quark model has been introduced fifty years ago mm as a realization 
of the SU{3) symmetry and it has been used with success for the description 
of many important properties of hadrons, as the existence of multiplets, their 
quantum numbers and the magnetic moments II 0). The idea of quarks as 
effective particles (Constituent Quarks) emerged very early [S] and was further 
developed with the introduction of the colour quantum numbers. 

Here we shall concentrate ourselves on Constituent Quark Models (CQM) 
for baryons. 

After the pioneering work of Isgur and Karl (IK) [H a series of CQM followed: 
the relativized Capstick-Isgur model (Cl) |7], the algebraic approach (BIL) |5], 
the hypercentral CQM (hCQM) [9l[Tn|[II|, the chiral Goldstone Boson Exchange 
model (yCQM) [121 [TS][TU[15] , the Bonn instanton model (BN) [THIITTIIT^fT^ 
and the interacting quark-diquark model [20] . 

All models reproduce the baryon spectrum, which is the first quantity to 
be approached when building a model for the baryon structure, but have been 
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widely used to describe baryon properties. In some cases the calculations re¬ 
ferred to as a CQM one are performed using a simple h.o. wave function for the 
internal quark motion either in a nonrelativistic (HO) or relativistic framework 
(rHO). 

The photocouplings for the excitation of the baryon resonances have been 
calculated in various models, among others we quote HO [Hj, IK [^, Cl [53], 
BIL [5], hCQM [53] (for a comparison among these and other previous ap¬ 
proaches see e.g. (SUES]). The calculations reproduce the overall trend, but 
the strength is systematically lower than the data. The fact that quite different 
models lead to similar results can be ascribed to their common SU (6) structure. 

As for the nucleon elastic form factors there are the calculations performed 
by BIL (SJ [53] with the algebraic method and by the Rome group |5S1 [53] 
within a light front approach based on the Cl model. The hCQM has been 
firstly applied in the nonrelativistic version with Lorentz boosts [solEI] and 
then it has been reformulated relativistically [351 [33] . A quite good description 
of the elastic form factors is achieved also using the GBE [Ml |3S] and the BN 
[33] models, both being fully relativistic. The same happens for the interacting 
quark-diquark model m, specially in its relativistic version m- 

A sensible test of both the energy and the short range properties of the 
quark structure is provided by the behaviour of the helicity amplitudes for 
the electromagnetic excitation to the baryon resonances. 

In the HO framework, there are various calculations of the transverse helicity 
amplitudes, among them we quote refs. [5iii55ii3HiEaiin], while a systematic 
rHO approach has been used by [41]. A light cone calculation, using the Cl |7] 
model, has been performed [53] and then successfully applied to the A [43] and 
Roper excitations [33]. For more recent light cone approaches, see ref. [33] and 
references therein. The algebraic method has been also used for the calculation 
of the transverse helicity amplitudes [5]. The hCQM, in its nonrelativistic ver¬ 
sion, has produced nice predictions for the transverse excitation of the negative 
parity resonances |45| and, recently, for both the transverse and longitudinal 
helicity amplitudes of all resonances having a sensible excitation strength [33] . 
The calculation of the helicity amplitudes in a relativistic hCQM is in progress 
and some preliminary results for the A resonance are now available [47] . He¬ 
licity amplitudes have been calculated also by the Bonn group, both for the 
nonstrange [331SH] and strange resonances [49] . 

The models have been applied also to the decays of baryons. The strong 
decays have been quite soon calculated with the IK model [55] and in its rel¬ 
ativized versions [sniiM]. There are also calculations in other models, namely 
BIL [55], GBE [53]. As for the hCQM, there are some preliminary calculations 
[53] . There are also calculations of the semileptonic decays of baryons in the 
BN model [55] . 

Finally we quote calculations of the axial nucleon form factors in the GBE 
[351153] and BN [33] . 
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2 A review of Constituent Quark Models 

2.1 Nonrelativistic approach 

The possibility of a nonrelativistic description of the internal quark dynamics 
was considered very early [5] after the introduction of the quark model. In this 
framework, one can introduce the three-quark wave function 1 * 35 , factorized 
according to the various degrees of freedom: 

^colour \spin ^flavour '^space- (I) 

In agreement with the Pauli principle, the wave function I'aq must be totally 
antisymmetric for the exchange of any quark pair. Baryons must be colour 
singlets and the corresponding wave function Ocoiour is by itself antisymmetric, 
therefore the remaining factors must be completely symmetric. Actually a sym¬ 
metric quark model has been formulated before the introduction of the colour 
quantum numbers and the symmetric three-quark states have been classified 

[571 [SH]. 

Early Lattice QCD calculations [SU] showed that the quark interaction can be 
split into a long range part, which is spin and flavour independent and contains 
confinement, and a short range spin-dependent one EO]. This means that one 
can assume the dominant part to be SU{6) invariant and the wave function of 
Eq. 0 becomes 

^'3, = ^colour ^SU(6) space • (2) 

In order to satisfy Pauli principle, the product 

^su(6) i’ space (3) 

must be symmetric and then both factors ^su{6) ijj space must have the same 

permutation symmetry, that is symmetric (S), antisymmetric (A) or one of the 
two mixed symmetry types (MS, MA), which are distinguished by the symmetry 
or the antisymmetry with respect to a quark pair. 

It should be reminded that each quark belongs to the fundamental SU{6) 
representation with dimension 6 and that with three quarks one can obtain the 
following SU ( 6 )-representations: 

S'I/( 6 ) : 6 (g) 6 (g) 6 = 20 0 70 0 70 © 56, (4) 

the corresponding symmetry type is, respectively. A, M, M, S. 

The spin and flavour content of each S'C/( 6 )-representation is well defined, 
since the three SU ( 6 ) representations can be decomposed according to the fol¬ 
lowing scheme 


20 

= U + ^8, 

(5) 

56 

= ^8 0 ^10, 

( 6 ) 

0 

II 

to 

+ ^8 + *8 + ^10. 

(7) 
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The suffixes in the r.h.s. denote the multiplicity 25'+1 of the 3q spin states and 
the underlined numbers are the dimensions of the 517(3) representations. This 
means for instance that the 56 representation contains a spin-1/2 5C/(3) octect 
and a spin-3/2 517(3) decuplet. 

The various baryon resonances can be grouped into SU ( 6 )-multiplets, the 
energy differences within each multiplet being at most of the order of 15% as 
in the case of 77 — A mass difference and of the splittings within the SU (3) 
multiplets. In Fig. we report the experimental non strange baryon spectrum, 
including only the three- and four- star states [6T] . The notation for the SU ( 6 )- 
multiplets is {d,L^), where d is the dimension of the 517(6)-representation, L 
is the total orbital angular momentum of the three-quark state describing the 
baryon and P the corresponding parity. An alternative but equivalent notation 
is Lf, where t is the symmetry type of the 517(6) representation. 

The fact that the 4— and 3—star non strange resonances can be arranged 
in 517(6) multiplets indicates that the quark dynamics has a dominant 517(6) 
invariant part accounting for the average multiplet energies, while the splittings 
within the multiplets are obtained by means of a 517(6) violating interaction, 
which can be spin and/or isospin dependent and can be treated as a perturba¬ 
tion. 


PDG 4* & 3* 



Figure 1: (Color online) The experimental spectrum of the non strange three- 
and four-star resonances m- The states are reported in columns with the same 
parity P and grouped into 517(6)-multiplets. 
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Table 1: Illustration of the features of various CQMs 


CQM 

Kin. Energy 

'^inv 


ref. 

Isgur-Karl 

nonrel. 

h.o. -I- shift 

OGE 

0 

Capstick-Isgur 

rel. 

string -kcoul-like 

OGE 

m 

U(7)B.I.L. 

IP 

vibr -|- L 

Giirsey-Rad 

0 

Hypercentral 

nonrel./rel. 

0(6): lin -|- hyp.coul 

OGE 

m 

Glozman-Riska 

rel. 

h.o. / linear 

GBE 

m 

Bonn 

rel. 

linear -|- 3 body 

instanton 

m 

quark-diquark 

nonrel./rel. 

linear -|- Coulomb 

spin-isospin 

m 


The various constituent quark models are quite different, but they have a 
simple general structure in common, since in any case, analogously to what 
stated above, the quark interaction V^q can be split into a spin-flavour inde¬ 
pendent part Vinv, which is S'{7(6)-invariant and contains the confinement in¬ 
teraction, and a S'f7(6)-dependent part 14/) which contains spin and eventually 
flavour dependent interactions 

V3q = V,nv + Vs f. ( 8 ) 

In Table we report a list of the Constituent Quark Models and their main 
features. The order is chronological. 


2.2 The Isgur-Karl model 

The Isgur-Karl model has some general features that are interesting also for 
other models, so it is worthwhile to devote to it some attention (more details 
can be found in [62]). 

The kinetic energy T is assumed to be nonrelativistic 

V = = Mtot-I-——- \-Tintr, (9) 

^ 2mj 2Mtot 


where is the total mass of the three quarks, P is their total momentum 
and the intrinsic kinetic energy is expressed in terms of the momenta px and 
Pp, which are conjugated to the Jacobi coordinates p and A, 


P= A= ^(fi-kr 2 -243) . 


( 10 ) 


In the case of non strange baryons, all quark have the same mass m and then 
T assumes the form 


( 11 ) 


in the case of quarks with different mass, the kinetic energy contains appropriate 
reduced masses El. 
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Table 2: The non strange three-quark states belonging to the h.o. shells 
up to N=2. The notation for the baryon resonances is X 212 J, where 
X=S,P,D,F,... denotes the pion wave in the decay channel, I and J are the 
isospin and the spin of the state, respectively. The numbers within parentheses 
are the masses of the 4- and 3- star resonances displayed in Fig. Q. The aster¬ 
isk in the first N=2 configuration reminds that the spin-isospin structure is the 
same as in the ground state but it has a radial excitation. 


N 

{d,Ln 

"8 

^8 

Ho 

no 

0 

(56,0+) 

Pii(939) 



7^33(1232) 

1 

(70,1-) 

5'ii(1535) 

S'ii(1650) 

^31(1620) 




iJi3(1520) 

7Ji3(1700) 

7^33(1700) 





7Ji5(1675) 



2 

(56,0+*) 

Pn(1440) 



7^33(1600) 

2 

(70,0+) 

Pii(1710) 

Pl3 

7^31 


2 

(56,2+) 

^13(1720) 



P3i(1910) 



Fi5(1680) 



7^33(1920) 






7^35(1905) 






7 ^ 37 ( 1950 ) 

2 

(70,2+) 

^3 

Pll 

P 33 




7^15 

Pl3 

P 35 





Fi5 






Fi7 



2 

(20,1+) 

Pll 






Pl3 





The confining interaction is assumed to be a harmonic oscillator (h.o.) 

= E ^ ^ - ^ 3 )^ (12) 

i<3 


which, in terms of the Jacobi coordinates, becomes 


Vno = \ k + A2), 


(13) 


with and = A^. The three-quark interaction is then given by two 

three-dimensional h.o. and the energy levels can be written as E = {3 + N)Tiuj, 
with N = 2n + Ip + l\, where n is a non negative integer number and Ip and l\ 
are the orbital angular momenta associated to the Jacobi coordinates. The h.o. 

parameter w is given by 

The general structure of the h.o. space wavefunction is 


V'ArLt(p,A) =C^P,v(p,A)e-i/2“^('’'+^') 1, flip) 11,(0^), 


(14) 
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where c? = 'JiKm/Ti, Cn is a normalization factor, Pff a polynomial of degree 
N and the spherical harmonics have to be combined to a definite total orbital 
angular momentum L; t (=A, M, S) is the symmetry type, the same as the 
SU{6) states. 

In Table we report the SU{6) states that can be assigned to the first three 
shells. All the states reported in Fig. Q fit very well into the scheme. The 
number of predicted states is however much larger than the observed 4- and 3- 
star states. The problem of such missing resonances is common to all CQMs and 
it has been suggested long time ago that some resonances may be observable 
in electroproduction experiments and not in strong interaction processes |22] . 
This statement is supported by the new states reported in the last edition of 
the PDG review [55] . 

In Tablej^ there is no room for the 3-star 1135(1930). In fact, the total spin 
5/2 can be obtained combining the L=1 total orbital angular momentum with 
the total spin 3/2 of the three quark, however the negative parity states belong 
to the 70-dimensional representation of SU(6), which cannot contain a A state 
with total spin 3/2. In order to describe this state and some new 2-star negative 
parity resonances [63] one should introduce the N=3 shell and the number of 
missing resonances will be highly increased. Similarly, the shells with N greater 
than 2 are necessary for the resonances with high spin values. 

An important observation regarding the h.o. spectrum is the level ordering, 
which, for any wo body potential is 0“*", 1“,0“*", while experimentally the 1“ 
states are in average almost degenerate with the first O’*' excitation. Moreover, 
the spacing between two shells is the same over the whole spectrum and the 
levels are highly degenerate, since the energy depends on the h.o. quantum 
number N only. 

In order to avoid the equal spacings and the degeneracy of the levels, in 
the Isgur-Karl model a shift potential U is added, which simply redefines the 
energies of the SU(6) states, without any attempt to diagonalize it. In this way 
the energies of the SU{6) configurations can be written as 


mt) 

— Eq, 

A(0+*) 

= Eq 2n — A, 

m~M) 

= Eo + 2n- A/2, 

mt) 

= Eo + 2n-2A/5, 

A(2+) 

= Eo + 2n- A/5 


= Eq -|- 20, 


where 


(15) 


n = ?ia; — ao /2 -I- 02/3, 


A =-5/4ao-b 5/302 - 1/304, (16) 


where the coefficients Om (ni=0,2,4) are determined by the moments of the U 
potential 


am = 3(^)3 / cPp{aprU{V2p)e-^^p\ 


( 17 ) 



No explicit form is assumed for the potential U, but the three coefficients am 
are used as free parameters to be fitted to the experimental spectrum. In this 
way also the position of the Roper N(1440) resonance is correctly described. 

Having assumed the space wave function as given by a h.o. three-quark 
potential, one can build the various SU{6) configurations to be identified, ac¬ 
cording to Table with the observed resonances. The states contained in each 
multiplet can be denoted as 

\B (18) 


where B = N, A for isospin 1/2, 3/2, respectively, S in the suffix is the total 
3q spin, J the spin of the baryon state, X=S,P,D,... according to the total 
3q orbital angular momentum and t is the symmetry type. For instance the 
nucleon is denoted by |7V ^£' 1 / 2 ) 5 , the Roper resonance is |A^ ^£*^ 2 ) 5 , where 
the asterisk means that the state is the first radial excitation of the nucleon, the 
A is IA ^£ 3 / 2)5 and so on. 

Since the quark interaction considered up to now is £C/( 6 ) invariant, the 
energies given by Eqs. (15) are common to all the states in any £1/(6) multiplet, 
at variance with the experimental spectrum (see Fig. [^. In order to describe 
the splittings within each multiplet, one has to introduce a £C/( 6 ) violating 
interaction 14/, which, in the case of the Isgur-Karl model is given by the 
hyperfine interaction, in the form proposed in ref. [60] 


Hhyp — 

i<j 


■ Si sm + 

'3 r3. 


-£*•£,], (19) 


where r) 


= n - ri. 


Eq. (191 is the spin dependent part of the One Gluon 


Exchange (OGE) interaction between two quarks, the spin independent part 
being a Coulomb-like term 1 /r^, which can be considered implicitly taken into 
account in the shift potential U. The structure of Eq. (191 is the same as the 
Breit-Fermi term in the higher order Coulomb potential for electrons in atoms. 
The OGE interaction is in principle valid for short interquark distances, however, 
it is used just for the determination of the form of the spin-dependent quark 
interaction and the strong coupling constant as is considered as a free parameter 
to be fitted to the N — A mass difference. 

The hyperfine interaction is diagonalized in the h.o. basis, using as unper¬ 
turbed energies the ones given by Eqs. (15). Its matrix elements, in the case 
of u and d quarks, are given in terms of the quantity (see also Appendix 2 of 

Ref. [55]) 


S = 


Aasa'^ 

3-\/27r; 


TTm^ 


( 20 ) 


which is substantially the N — A mass difference and can be fixed to about 
300 MeV. As for the remaining free parameters, the h.o. constant is fitted to 
the proton r.m.s. radius, obtaining = 1.23 fm^, the other parameters are 
determined by comparison of the theoretical spectrum with the experimental 
one. The resulting description of the spectrum is quite good, both for non 
strange and strange resonances | 6 ]. 
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An important consequence of the introduction of the hyperfine interaction is 
that the baryon states are superpositions of SU{6) configurations. For instance, 
the nucleon is expanded as 


1-^) — 0,s\N‘^ S 1 / 2 ) s + ^ 1 / 2)3 + S 1 / 2 ) M + 0,d\N'^Di/2)m, (21) 


with as = 0 . 931,05 = —0.274,om = —0.233,od = —0.067 (55]; the asterisk in 
the second term of Eq. (21) means that the spin-isospin part is the same as the 


first term, but the space part corresponds to a radially excited wave function. 
The Roper resonance has a similar expansion, with the dominant component 

given by \N‘^S[^^)s' Os = 0.281, a( = —0.960, om =- 0.003, an = —0.001. 

The A resonance is given by 


|A) = 65|A453/2)s + 6 's|A"^ 3 */ 2 )s + bnlA%^ 2 }s + b'^IA^n3/2}M, (22) 

with 65 = 0 . 963,65 = 0 . 231,651 = —0.119, 6 ' 5 , = 0.075. It is well known that 
with pure SU(6) configurations the E2 electromagnetic N — A transition van¬ 
ishes |64| . However, because of the hyperfine interaction, the A state acquires 
a non-zero D-wave component and then a small quadrupole strength arises 
[67l l 66 ] . The theoretical estimate of the ratio 


R = - 


Ge2 

Gmi 


(23) 


is about —0.02 [661l67j . which compares favourably with the experimental value 

E3]- 

The behaviour of the nucleon electromagnetic form factors in the Isgur-Karl 

_ 

model is dominated by the Gauss factor e ^ therefore it is too strongly 
damped for medium-high values of the square momentum of the virtual photon 
. On the contrary, the neutron charge form factor is nicely described and this 
is due to the hyperfine interaction [68) . In fact, if the nucleon state is the sym¬ 
metric SU{6) configuration \N^Si/ 2 )sj the charge form factor is proportional to 
the total charge of the three-quark system; the hyperfine interaction introduces 
in the nucleon state a mixed symmetry component \N‘^Si/ 2 )m, giving rise to a 
non zero charge form factor for the neutron. 


2.3 The Capstick-Isgur model 

This model [7] is the extension to the baryon sector of the relativized model for 
mesons formulated in ref. [69) . The three-quark hamiltonian is written as 

H = T + V3g, (24) 

where T is the relativistic kinetic energy 

T = Y. + (25) 
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the three quark potential Vsq is separated into two terms, according to Eq. 

In the nonrelativistic limit 

V3q ^ E,, + Vsd. (26) 

where the spin dependent interaction is 

Esi — ^string + Vcoul- (27) 


The first term is the three-body adiabatic potential generated by the quantum 
ground state in a F—shaped string configuration and provides confinement; 
Vstnng is given by [7] 


3 

^string — ^ ^ ( \^i '^junction\ 

i=l 


(28) 


where Cqqq is an overall constant energy shift and b is the string tension. For 
practical purposes, Vstring is split into two-and three-body effective terms 




string 


= Ca 


i<j 




(29) 


where 


2 

^3b — ^ junctionl f ^ ^ )' 

i<j 


(30) 


The parameter f is chosen to be 0.5493 m in order to minimize the expectation 
value of Vsb in the h.o. ground state of the baryon. In this way Vsb is a small 
correction and can be treated perturbatively and Vstring becomes very close to 
1/25 Tij. 

The potential Vcouh in the nonrelativistic limit, is given by 


Vcoul — ^ ( 


i<j 


basing) ^ 

Q^. . ’ 

Ol ij 


(31) 


in ref. [7] the momentum (or space) dependence of the strong coupling constant 
oisifij) is properly taken into account. 

The spin dependent potential Vsd-, again in the nonrelativistic limit, is 


Vsd — Vhyp -f Vso: 


(32) 


where Vhyp is the hyperfine interaction of Eq. (191 and Vso is a spin-orbit in¬ 
teraction containing in particular a Thomas precession term. Please note that 
the sum of Vcoul and Vhyp, together with the corresponding Thomas precession 
spin-orbit, derive from the nonrelativistic limit of the OGE interaction. 

In order to avoid the nonrelativistic approximation one has, according to 

the discussion reported in ref. [7], to introduce appropriate factors and to 
smear the interactions over a two quark distribution 


Pij(fi -rj) = 


r3/2 




(33) 
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in particular for the contact term in the hyperfine interaction. In this way the 
factors rrii are substituted with the corresponding energies and the momentum 
dependence of the interaction is taken into account. 

The three-body equation, with the relativistic kinetic energy, is solved by 
means of a variational approach in a large h.o. basis. The result is a good 
description of the baryon spectrum, including strange, charm and bottom reso¬ 
nances [7]. 

As already mentioned in the introduction, the model by Capstick-Isgur has 
been successfully applied to the calculation of the electromagnetic amplitudes 
for the transitions to the A [32] and Roper resonances [43] . 


2.4 The U(7) model 

The typical feature of this model |8] is to describe the state of a three quark 
system by means of a group theoretic approach. 

In order to describe the space degrees of freedom, the model uses the method 
of bosonic quantization, similarly to what has been done in the Interacting Boson 
model in nuclear physics mi and in molecular physics m as well. The idea 
is to consider a string-like model with a Y-shaped configuration, in which the 
vectors fi{i = 1,2,3) denote the end points of the string configuration. To this 
end, one introduce s tw o vector boson operators defined in terms of the Jacobi 
coordinates of Eq. (101 p, A, together with their conjugate momenta p'p, p\ 


bp^m — 


T iVpjm) 1 


hi. 





(Ayn T fPA,m), 


t 


1 

71 


{\m ‘^PX,m') 


(34) 

(35) 


(with m = —1,0,-fl) and an auxiliary scalar boson s, sb The bilinear forms 
Gaa' = cica', where (a = 1,...7) is one of the seven creation operators, 
generate the Lie algebra of U(7). The choice of U(7) is in agreement with the 
usual prescription that any problem with ix space degrees of freedom should be 
written in terms of the Lie algebra U(iy + 1) [73] and all the states are assigned 
to the totally symmetric representation [N] oi U{i'+ 1), N being the maximum 
number of shells. The physical states can be constructed by applying a suitable 
product of boson operators to the vacuum state 


(36) 

AC is a normalization factor. 

The squared mass operator is expressed as the most general combination 
of Gaa'i with the condition of being at most quadratic, preserving angular 
momentum and parity and transforming as a scalar under the permutation 
group. The general form of contains several models of baryons structure, 
including single particle (i.e. h.o.) and collective string models. The calculations 
are performed choosing the latter model,with the consequence that contains 
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also a term of the type ss + h.c., which causes a spread of the wave function 
over many h.o. shells. In order to make more transparent the interpretation of 
the results, the mass operator is rewritten in terms of vibrational and rotational 
contributions to the baryon spectrum, using a procedure already introduced for 
the Interacting Boson Model m 

(37) 

in this way the baryon excitation spectrum is determined by vibrations and 
rotations of the string-like configuration. There seems to be no evidence of 
excitations due to the term therefore it is omitted and the remaining 

two terms, according to the discussion reported in [S], are simplified obtaining 
the mass formula 


— AIq N\^KiTlu “b “b CtT, 


(38) 


where /ti,K 2 ,a are free parameters, L is the total orbital angular momentum 
and n„, n«, Uu, are the eigenvalues of number operators of the type c^c, labeling 
the vibrational energy levels. 


Eq. (38) describes the SU{6) invariant part of the interaction. The splittings 


within the multiples are introduced with reference to the internal part of the 
state (see Eq. 0). The corresponding algebraic structure is 

Gi=SUc{3)<8)SUs{2)(S)SUf{3), (39) 

describing the colour, spin and flavour degrees of freedom, respectively. Baryons 
are colour singlets and then only the spin-flavour degrees of freedom contribute 
to the energy splittings. The mass squared operator is written in a Giirsey- 
Radicati form m 

a[C2{SU{%)sf) - 45] + b[C2{SU{i )/) - 9] + 


= 


+b'[C2{SUi{2)) - |] + 6"[Ci(t/y(l)) - 1] + 


(40) 


+b"'[C2{UY{l)) - 1] + c[C2{SUs(2)) - fj. 

The quantities denoted as Cn{X),n = 1,2 are the Casimir operators of the 


Group X; the constants in Eq. (40) are chosen in order that each term vanishes 
in the nucleon ground state. Eor non strange baryons the hypercharge Y is 
equal to 1 and the b and b’ terms can be grouped into a single term, therefore 
one can use the simplified form 


= a[C2iSU{6),f) - 45] 


b[C2iSU{3)f) - 9] + c[C2{SUsi2)) - -]. 


(41) 


The seven parameters in Eqs. (38) and (41) are obtained fitting the non 


strange baryon spectrum and the results are very good. In the model the num¬ 
ber of shell is not limited, therefore one can describe well also the four star res¬ 
onances with higher values of the spin, such as G19,7719, G19,/111 and 7/311. 
The extension of the model to the strange baryons is presented in ref. m- 
The model allows to calculate also the elastic nucleon form factors [55] and 
the electromagnetic transition amplitudes for photo-and electroproduction jS] 
[26] , provided that a form for the charge distribution along the string is assumed. 
In this way both the elastic and inelastic form factors are adequately described. 
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2.5 The Goldstone Boson Exchange Model 

The model is based on the consideration that QCD exhibits an approximate 
chiral symmetry which is spontaneously broken |12| . As a consequence of such 
spontaneous symmetry breaking, quarks acquire an effective mass and Gold- 
stone bosons emerge, which are indentified with the pseudoscalar meson octet. 
Therefore it is assumed that baryons are considered as a system of three con¬ 
stituent quarks with an effective quark-quark interaction Va^, which is split into 
parts according to Eq. While different forms are assumed for Vmv in the 
various versions of the model, the spin-flavour part is always chosen as an ex¬ 
change of pseudoscalar (Goldstone) bosons between two quarks. The simplest 
form of this chiral interaction can be written as [12] 


i?v 


i<j 


V{u,)xl 




(42) 


where \f are the flavour Gell-Man matrices and ai the quark spin operators. 
The interaction has a Yukawa form containing a spin-spin and a tensor part. 
The spin-spin part is given by 




47r 3 4:7711771 j 


1 ^ 

Si ■ SjXf ■ X^[^JS -47r(5(rij)], 


(43) 


where rrii are the quark masses and P labels the exchanged boson of mass /x; the 
5 is actually smeared out by the finite size of quarks and mesons. The flavour 
structure of the quark-quark interaction is then 


Voctetinj) = 'yi V^{r,j)Xi ■ X' 


a—1 


7 

6=4 


V^(r.,)A^A5 + E,(r.,)A^A«. (44) 


In ref. [T^] Vinv is assumed to be a h.o. potential and the boson exchange 
interaction is considered in the chiral limit, in which case the masses of all the 
three quarks are equal and 14- = Vk = Ej. Treating the interaction as a pertur¬ 
bation, the masses of the baryons can be expressed in terms of a limited number 
of radial integrals, which are used as parameters in order to fit the experimen¬ 
tal values. Already in this simplified approach, a reasonable description of the 
spectrum is obtained, in particular the boson exchange interaction leads to the 
correct ordering between the excited O’*" and the first negative parity levels. As 
a further improvement also deviations from the chiral limit and the contribution 
of the tensor part of the boson exchange interaction are considered. 

A substantial improvement of the model has been started in ref. |13j , in the 
sense that the confinement interaction is assumed to be linear and V^q is inserted 
into a Faddeev equation for the three quark system to be solved numerically. 


The Vsf is given by Eq. (44), to which the exchange of the singlet meson (rj') 
is added. In the first application, only nucleon and A states are considered 
[13 m]. The interquark potential is then: 


— ^octeti'^ij) T ^singlet Ct. 


ij) 


(45) 
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( 46 ) 


in Voctet only the tt and 77 potentials contribute and the singlet potential is 

2 ^ 


The quark-eta coupling constant is assumed equal to the quark-pion one, de¬ 
duced from the pion-nucleon interaction. Keeping the meson masses equal to 
their physical values, there remain only four free parameters, namely the two 
parameters which determine the smearing of the <5 function in Eq. (431, the rj' 
coupling constant and the strength of the linear confinement C. The result is 
a good description of the 14 lowest N and A states, respecting the ordering 
displayed by the experimental spectra. 

A unified description of both non strange and strange baryons is finally 
achieved in m, where the interaction of Eq. is used together with a rela¬ 
tivistic kinetic energy as in Eq. (25). The three-quark wave equation is solved 
by means of a variational approach. Again a quite satisfactory description of 
the low-lying light and strange baryons is achieved, respecting in particular the 
already mentioned relative ordering of the positive and negative parity states. 


2.6 The Bonn model 


The authors start from the consideration that the nonrelativistic approach seems 
to be completely inadequate for the description of the internal motion of quarks 
with small constituent masses and therefore they introduce a relativistic formu¬ 
lation [TB] . 

The relativistic formulation is performed within quantum field theory and is 
based on the six-point Green’s function Gx describing three interacting quarks. 
The infinite series of Feynman diagrams, necessary in order to describe a bound 
state, is rearranged in the same way used by Bethe-Salpeter for the two-particle 
case. The result is that Gx obeys to an integral equation containing two irre¬ 
ducible kernels and describing the two- and three-particle interac¬ 

tions, respectively. Introducing the momentum space representation Gp of Gx-, 
the integral equation can be written in concise form as 

Gp = Gqp — iGopKpGp, (47) 


where Gqp is the three quark propagator and Kp the total kernel, or equivalently 
as 

(G-^ + iKp)Gp = I, (48) 

showing that Gp is the resolvent of the pseudohamiltonian Hp 

Hp = {G-^+tKp). (49) 


The idea is to extract from Gp the baryon contributions, meant as real 
bound states of three quarks with positive energy _j_ ]\/[2^ 
is expanded in a Laurent series, which, near a pole, gives 


Gp{pp,P\]p'p,p'x) 


. Xp{Pp,Px)Xp{Pp,Px) 

* P^-M^+ie 


(50) 
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The quantity XpiPp^Px) is the Fourier transform of the Bethe-Salpeter (BS) am¬ 
plitude Xp(a;i, a; 2 , cca) for the bound state |P), defined as transition amplitudes 
between the state |P) and the vacuum |0) 


Xp{xi,X2,Xz) = a2ix2)'^ a3{x3))\P) , 


(51) 


where is the quark held and Oi denotes the Dirac, flavour and colour 

indices. Thanks to translational invariance, one gets a BS amplitude Xpip, A) 
which depends on relative coordinates only and, in momentum space, on the 
conjugate momenta Pp,p\\ here p and A are tetravectors, whose spatial part 
coincide practically with the Jacobi coordinates dehned in Eq. (10). 


The factorization property of the pole residue allows then to extract the BS 
amplitude Xp{PpiPx)^ which satishes the equation 


Xp Apj 


(52) 


where has, as mentioned before, two- and three-body contributions 
and respectively. 


In principle the BS equation (52) allows a covariant description of baryons 


as bound states of three quarks in the framework of QCD. However it cannot 
be used practically because the single quark propagators and the kernels K^'^\ 
are only formally dehned in perturbation theory as an inhnite sum of Feyn¬ 
man diagrams and are not deducible from QCD. Moreover, the dependence on 
the relative energy (or relative time) leads to a complicate analytical structure. 
Therefore an adequate parametrization is necessary, to be introduced after hav¬ 
ing having obtained a six-dimensional reduction of the full eight dimensional BS 
equation, the so called Salpeter equation, trying to preserve the covariance of 
the theory and keep it as close as possible to the quite successful nonrelativistic 
quark model. 

This can be achieved following the lines of what has already performed in the 
covariant quark model for meson case using an instantaneous qq BS equation 
EH- To this end, the quark propagators are assumed to be given by their free 
forms with effective constituent masses. Moreover the kernels and are 
approximated by effective interactions which are instantaneous in the baryon 
rest frame, that is 


Kp\pp,P\\p’p-,p'\) = yP{Pp,Px',P'p,p'x), (53) 

Mp+p^Pp./p) = yp\PpX,^p), (54) 

where P = {M,0). This instantaneous approximation can be formulated in 
a covariant way following the method proposed in ref. [78) . The reduction to 
the six dimensional Salpeter equation can be performed more easily if only the 
three-body kernel is present; the introduction of the two body kernel is possible 
provided that it is substituted by a suitable effective interaction 14//. In any 
case the reduction to the Salpeter amplitude is achieved by means of an 
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integration over the energy variables, a procedure which, thanks to Eqs. (531 
and (54) affects the BS amplitude only 


<^m{Pp,P\) = J ^^Xm{p%pI,Pp,P\), (55) 

where xm is the BS amplitude in the rest frame. 

Finally the Salpeter equation is written in a Hamiltonian formulation 

UmK = (56) 

where A is a projector operator over positive energy states. 

In order to perform explicit calculations of the baryon spectrum one has to 
assume some specific form of the hamilton operator Hm ■ In agreement with the 
previous discussion, Hm contains two-and three-body potentials m- 

The confining three-body potential is chosen within a string-like picture, 
where the quarks are connected by gluonic strings (flux tubes) and the potential 
increases linearly with a collective radius r^q 

^2, ra) ~ r3q{fi,f2,r3)- (57) 

There are three different ways to define r^q [T7]. The first one is the T-type [7S] 

3 

rsq = ry = min'^\f, - fo\, (58) 

fo is the position where the fiux tubes can merge and is chosen in order to 
minimize r 3 g. A second possibility is given by the so called A-type 

3 

C 39 = rA =y]|ri-ril; (59) 

i<j 

rescaling ta by a factor f one gets a good approximation of ry m and ref¬ 
erences quoted therein), provided that 1/2 </ < l/-\/3. The third choice is 
provided by the hypercentral one [5] 

rsq = rH = a/A+T. (60) 

Having tested that the structure of spectra depends only slightly on the choice 
of r 3 q, the authors of ref. [T7j use the three-body A-shape string potential rising 
linearly with ta, which provides the SU{6) invariant part of the three-quark 
interaction. 

The two-body potential is taken from the instanton interaction introduced 
by 't Hooft [HD] in the SU{2) case (extended to SU{3) in ref. [HT]), in which the 
S term is smeared with an effective range A. Such an interaction acts only on 
flavour antisymmetric states and therefore it does non act on A states, thereby 
leading to a iV — A mass splitting. 
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The model hamiltonian depends on seven free parameters, which are used to 
describe both the non-strange m and strange [H] baryon sector. The results 
are quite satisfactory. 

The model has been successfully applied to the description of the elastic 
form factors [3S1135] , the helicity amplitudes for both the nonstrange [15] and 
strange resonances [IHI, the semileptonic decays of baryons jSSj and the axial 
form factors [55] . 


2.7 The interacting quark-diquark model 

The Interacting quark quark model [20] and its relativistic version [52] give 
a good reproduction of the spectrum, moreover they have much less missing 
resonances than a normal three quark model. In particular, we report here the 
rest frame mass operator of the Relativistic quark-diquark model : 


M = Eo + 



(61) 


where Eq is a constant, Mdir{r) and respectively, are the direct and the 

exchange diquark-quark interaction, mi and m 2 stand for the diquark and quark 
masses, where m-i is either ms or mAv according if the mass operator acts on 
a scalar or an axial vector diquark, and Mcoatij") is a contact interaction. The 
direct term is a Coulomb-like interaction with a cutoff plus a linear confinement 
term 

Mdn = -^(l-e-^^)+/3r. (62) 

A simple mechanism that generates a Coulomb-like interaction is the one-gluon 
exchange. One needs also an exchange interaction. This is indeed the crucial 
ingredient of a quark-diquark description of baryons and has the form 

Alex(r) = (-l)'+^2Ae“‘^’' [As(si • S 2 ) + Ajiti ■ h) + As/(si • S2)(ri • ^" 2 )] , (63) 


where s and t are the spin and the isospin operators. Moreover, we consider a 
contact interaction 


Mcont(r) = 


/ WlW2 \ 

V E1E2 j 


1 / 2 +e 


7r3/2 


2 ^ .c 

\ E 1 E 2 


(64) 


where Ei = q"^ + m.|(z = 1,2), 6,77 and D are parameters of the model. The 

relativistic Interacting quark-diquark model is a relativistic version of the In¬ 
teracting model of Ref. [20] . The Interacting quark-diquark model hamiltonian 
is 

H = Eq -\- ^ - \- Pr + [B + 0(5^ o]<5si,i 

+ (—l)*^^2Ae [(si • S 2 ) -l- {ti ■ t2) + (si • S 2 )(ti • O)] , (65) 

where si and 5*2 are the spin of the quark and of the diquark respectively, while 
ti and <2 the the same for the isospin. The contact interaction C5l,o acts only 
on the spatial ground state, while the ^si.i on the axial diquark. 
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Figure 2: (Color online) The experimental spectrum of the non strange three- 
and four-star resonances m in comparison with the results of the interacting 
quark-diquark model |20) . 


3 The hypercentral Constituent Quark Model 


3.1 The hyperspherical coordinates 


The starting point of the hypercentral Constituent Quark Model (hCQM) is the 
introduction of the hyperspherical coordinates |HSl HH Hi] , which are given by 
the angles tip = {6p,(j)p) and together with the hyperradius, x, 

and the hyperangle, defined in terms of the absolute values p and A of the 
Jacobi coordinates of Eq. (10) 


X = 


^ = arctan ^ • 
A 


( 66 ) 


The hyperradius x is a collective variable, which gives a measure of the dimension 
of the three-quark system, while the hyperangle ^ reflects its deformation. 

Using these variables, the nonrelativistic kinetic energy operator of Eq. 
after having separated the c.m. motion, can be written as 


AN 9 ^ 5 9 

2m^ 2m^dx^^xdx^ x^ ^ 


(67) 
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The grand angular operator L^(n) = L'^{ilp,nx,^) is the six-dimensional gen¬ 
eralization of the squared angular momentum operator and is a representa¬ 
tion of the quadratic Casimir operator of the rotation group in six dimensions 
0(6). Its eigenfunctions are the so called hyperspherical harmonics (h.h.) [55] 


the grand angular quantum number 7 is given hy j = 2n + Ip + lx, where n is 
a nonnegative integer and Ip, lx are the angular momenta corresponding to the 


Jacobi coordinates of Eq. (10). 


The h.h. describe the angular and hyperangular part of the three-quark wave 
function and are written as 1851 


= F ™ flip) F,,„,(0;,) (^). 


(69) 


where the hyperangular functions (^) are given in terms of trigonometric 

functions and Jacobi polynomials [55] . 

The h.h. form a complete orthogonal basis in the space of the functions of 
Up, fix, ^ and then any three-quark wave function can be expanded as a series 
of h.h. 


^/)^(a:) (70) 

the hyperradial wave function 'tp^{x), depending on the hyperradius x only, is 
completely symmetric for the exchange of the quark coordinates. 


3.2 The hypercentral approximation 


An expansion similar to Eq. (70) is valid for any quark interaction, the Hrst 
term depending on the hyperradius x only: 


V{n,) = V{x) + • • •. (71) 

Retaining the first term only one gets the so called hypercentral approxima¬ 
tion. Such approximation has been applied with success to the description of 
few-nucleon systems | 86 l [87] , while in the baryon case it has been shown that 
the matrix elements of the currently used two-body qq potentials in the 3q space 
exhibit an almost perfectly hypercentral behaviour [ 88 ] . 

In the hypercentral approximation, the three-quark potential depends on 
the hyperradius x only and therefore it has a three-body character, since the 
dependence on the single pair coordinates cannot be disentangled from the third 
one. The possibility of three-body forces is strictly related to the existence of 
a direct gluon-gluon interaction, which is one of the fundamental features of 
QCD. The diagram shown in Fig. [^a is the lowest order one leading to a non 
vanishing three-body interaction among quarks in a baryon, but of course many 
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t * N 



Figure 3: a) The lowest order diagram leading to a non zero three-body in¬ 
teraction of quarks in a baryon. b) The -Y-shaped string configuration, c) The 
A-shaped string configuration. 


others can be considered. A three-quark mechanism is considered also in flux 
tube models, which have been proposed as a QCD-based description of quark 
interactions [89], leading to the Y-shaped three-quark configuration of Fig. Ep 
izsi, besides the standard A-like two-body one of Fig.[^. Furthermore, a Born- 
Oppenheimer treatment of the confinement potential in a QCD motivated bag 
model leads quite naturally to three-body forces (55] [SD] which increases linearly 
with some ’collective’ radius. 

In the hCQM the three-quark interaction is assumed to be hypercentral 

F 3 ,(p,A) = V{x), (72) 


as a consequence, the three-quark wave function is factorized 


V’3(3(PjA) — (flp, ^); 


(73) 


the hyperradial wave function is labeled by the grand angular quan¬ 

tum number 7 defined above and by the number of nodes v. The angular- 
hyperangular part of the 3q-state is completely described by the h.h. and is the 
same for any hypercentral potential. The dynamics is contained in the hyper¬ 
radial wave function which, because of the factorization Eq. (73), is 

obtained as a solution of the hyperradial equation 


(f 5 
dx'^ 


d 

X dx 


7(7 + 4 ) 


X‘^ 


-] = - 2m [E - V 3 q{x)] 'tp^„{x). (74) 
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The Eq. (74) can be solved analytically in two cases. The first one is the 
six-dimensional harmonic oscillator (h.o.) 


^ i fc (rl - fjf = = Vh.o{x), (75) 

i<j 


which is exactly hypercentral. The eigenvalues are given, as already mentioned, 
by £1 = (3 -I- where N can be written as N = 2]^ + j and the hyperradial 

wave functions are reported in Appendix A. 

The second analytical case is given by the hyperCoulomb (hC) potential 

[iniiniiisiiMiini] 

Vhye{x) = -^. (76) 

The eigenvalues of the hyperCoulomb problem can be obtained by generaliz¬ 
ing to six dimensions the calculations performed in three dimensions, obtaining 


2 

T m 


(77) 


where n = TV -I- | is the principal quantum number and N = 7 - 1 - 1 /, where 
V = 0 , 1 , 2 , ... is the radial quantum number that counts the number of nodes 
of the wave function. 

The fundamental reason why the three-body problem with h.o. or hC inter¬ 
action is exactly solvable is that they have a dynamic symmetry, 17(6) and 0(7) 
respectively. 

The dynamic symmetry 0(7) of the hC problem can be used to obtain the 
eigenvalues using purely algebraic methods. The hyperCoulomb Hamiltonian 
can be rewritten as [32] 


H 


2 

T m 

2 [02(0(7))+ f] ’ 


(78) 


where 02 ( 0 ( 7 )) is the quadratic Casimir operator of 0(7). It can be shown [Si] 
that the eigenvalues of 02 ( 0 ( 7 )) are given by {v + ^){v -|- 7 -|- 5), obtaining 


E 


2 

T m 

2 [( 1 / -I- 7 )(i/ -I- 7 -b 5 ) -b 25/4] ’ 


(79) 


which coincides with Eq. ( 0 . 

The eigenfunctions of Eq. ( |74| with the hyperCoulomb potential can be 
obtained analytically and are |M) 




(2g)" 

(27 + 2v + 5){v -b 27 -b 4 )!^ 


i2gxre-^^ Ll'^+\2gx), (80) 


where for the associated Laguerre polynomials the notation of Ref. |3S| is used 
and g = . The explicit expression of the hyperradial wave functions are 

obtained in ref. |94| and reported in Appendix A. 
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Figure 4: (Color online) Qualitative structure of theoretical spectra for the h.o. 
(left) and for the hC potentials (right), up to the first three shells. The energy 
units are arbitrary and different for the two potentials. 


A complete solution of the hyperCoulomb problem, using the SO(7,2) dy¬ 
namical group, has been worked out in ref. |96| . where the states, the elastic 
and inelastic form factors have been also developed. 

The hC potential has important features miiaiMi- The energy eigenvalues 
depend on i/ -b 7 and then the negative parity states are exactly degenerate 
with the positive parity excitations, as is shown in Fig. The observed 

Roper resonance is somewhat lower with respect to the negative parity baryon 
resonance, at variance with the prediction of any S'[/( 6 )—invariant two-body 
potential, therefore the hC potential provides a better starting point for the 
description of the spectrum. The spectrum of Fig. Q shows that within the 
first three shell the hC potential exhibits two more states with respect to the 
h.o. The first extra level has positive parity and contains a further N and A 
state, thus enhancing slightly the number of theoretical states, but the second 
one has negative parity and allows to insert the recently observed states, already 
mentioned in Sec. 2.2. 

Another interesting property of the hC potential is that the form factors 
calculated with its wave functions have a power-law behaviour uni HI HD, 
leading to an improvement with respect to the widely used harmonic oscillator, 
for which the form factors are too strong damped for increasing momentum 
transfer. 

3.3 The hypercentral Constituent Quark Model 

The hC potential has interesting features, however it is not confining. In the 
hCQM, the SU{6) conserving part of the potential is then assumed to be the 
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Figure 5: (Color online) Qualitative structure of theoretical spectra for the 
potential of Eq. (14), up to the first three shells. The energy units are arbitrary. 


sum of the hC interaction and a linear confinement term ISE] 


hint; — T 

X 


(81) 


the hyperradial equation (74) must be solved numerically and the presence of 


the confinement removes the degeneracies typical of the hC potential, as shown 


in Fig. (141, however the general structure is only slightly modified with respect 


to the hC potential. 


The structure of the potential of Eq. (81) is formally similar to the quark- 


antiquark Cornell potential m widely used for the description of mesons. It 
is noteworthy that Lattice QCD calculations [98] are able to reproduce the 
Cornell potential, as it is seen in Fig. ([^, where the LQCD results for static 
quark-antiquark pairs in the SU{S) limit are reported. The model potential 
of Eq. (81) can then be considered as the hypercentral approximation of the 


Cornell potential. 

Having chosen the form for the hypercentral potential, the solutions of the 


hypercentral equation (74) produce a series of wave functions ip^,y{x) and then 


one can build up the model SU{6) states. Taking advantage of the fact that 
also the h.o. potential is hypercentral, one can start from the states of the 
Isgur-Karl model jS], express them in terms of the hyperspherical coordinates 
and substitute the h.o. hyperradial wave functions with those determined by 


the potential of Eq. (81). The complete SU{6) configurations for non strange 


baryons are reported in Appendix B. 
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Figure 6: (Color online) The quark-antiquark potential calculated in LQCD 
[98] for static quarks in the SU{3) limit. The constant /3 is the inverse QCD 
coupling and tq ~ 0.5 fm. 


In order to complete the model hamiltonian one has to add a term violating 
the SU{6) symmetry. In hCQM such term is chosen to be the standard hyperfine 
interaction (SOI IS] of Eq. (U^. The hamiltonian for the three quark system in 
the hCQM is then 


HhCQM — 3to -|- ~ X Hfiyp. (82) 

In this way the baryon states are superpositions of the SU{6) configurations 
reported in Appendix B. 

Using the notation introduced with the Isgur-Karl model, the nucleon state 
can be written as 


\N)= as\N^Si/2)s + a's\N^Sl/^)s 




+aM|A^^>S'i/2)M + az>|A^^T*i/2)M) 


(83) 


with as = 0.976, Og = —0.196, Og = —0.043, om = —0.051, od = —0.070 [9l [TT]: 
the asterisks in the second and third term of Eq. (831 mean that the spin- 


isospin part are the same as in the first term, but the space part corresponds to 
hyperradially excited wave functions. One should not forget that in the hCQM 
there are two hyperradial excitations of the nucleon within the first three shells 
(see Fig. (|^. The Roper resonance has a similar expansion, with the dominant 
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Figure 7: (Color online) The spectrum obtained using the hCQM hamiltonian 
of Eq. (821. The free parameters are fitted to the experimental values of the 4* 
and 3* resonances reported in the PDG m- 


component given by “s = 0.183, a), = 0.967, a" = 0.185, om = 

0.021, a_D = —0.176. The A resonance is given by 

1^) = I ^"*’< 53 / 2)5 + b'g\A'^Sy2) S + ^r>l^^'S'3/2)s + ^dI^^^3/2)m, (84) 

with bs = 0.976, b'g = 0.166,6^ = -0.042, bo = 0.026, b'jy = -0.132. 

The quark mass is taken to be 1/3 of the proton mass, as prescribed in 
order to reproduce the proton magnetic moment (see e.g. |62jb In this way in 
the model there are only three free parameters: r, a and the strength of the 
hyperfine interaction, the latter being mainly determined form the N — A mass 
difference. The parameters are found by htting the energies of the 4* and 3* 
non strange baryons relative to the nucleon and are given by 

a=1.16/m-2, t = 4.59. (85) 

The resulting spectrum is reported in Fig. 0- 

Having fixed the parameters of the three-quark hamiltonian from the spec¬ 
trum, the baryon states are completely determined and can be used for the 
calculations of various properties. In the rest of the paper, the results obtained 
by the present nonrelativistic hCQM are given by parameter free calculations, 
that is they are predictions. 


3.4 An analytical model 


The hyperradial equation with the hypercentral potential of Eq. (81) cannot 
be solved analytically unless the linear confinement term is treated as a per¬ 
turbation [93l [94] . This situation is reasonably valid for the lower states since 
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Figure 8 : (Color online) The spectrum obtained using the analytical model. 
The free parameters are htted to the experiment! values of the 4* and 3* reso¬ 
nances reported in the PDG m- 


they are confined in the low x region where the hC term is dominant. With 
this assumption, the perturbative contributions to the energies of the states are 
determined by the integral 


which is given by [Ml [M] 

1 i/\ 

2g {iy + 2-f + 5){iy + 2-f + 4)!^ L- + V- “T - / “T ^ 2)12 ’ 

(87) 

where 

V — 2 < a < V. ( 88 ) 

The energy eigenvalues are then 


dx \'ip-yv{x)'\^ 


( 86 ) 


ro T/' I o r( 7 -l -27 + CT) 


where, according to the definition introduced in Sec. 3.2, n is given by z^-|- 7 -|- 5 / 2 . 
This formula, to which a constant Eq should be added, is used to reproduce the 
energies averaged over the states with the same quantum numbers v, 7 . In the 
fitting procedure, the Roper resonance is not taken into account because the 
presence of the confinement pushes it upwards with respect to the negative 
parity resonances. Assuming a quark mass about 1/3 of the proton mass, the 
result of the fit leads to [Ml El] 

Eo = 2.152 GeV, r = 6.39, a = 0.148 fm■^ (90) 
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the latter two values are not much different from the ones of the non perturbative 


analysis reported in Eq. (85). 


In order to describe the splittings within the SU{6) multiplets, one has to 
add a SU{6) violating interaction to the potential of Eq. (81). Eor simplicity, 
such interaction is assumed ESI 


to contain a spin-spin term 


V^{x) = A 


E 

z<j 


= A 


[2 


where S is the total spin of the 3-quark system, and a tensor interaction 


B 4 E 


Kj 


(a, ■ (rl - rj)) (d) ■ (rl - f})) 1 ^^ ^ 

-— - O 


n - rP 


(91) 


(92) 


The parameters of the spin-spin interaction Eq. are determined from 
the N{939) — A(1232) and the fV(1535) — A(1650) splittings, obtaining 
A = 140.7 MeV and /3 = 1.53 fm“^. We observe that, at variance with the 
OGE hyperfine interaction, the spin-spin term has a non-zero range, a feature 
that seems to be necessary since the hC wave functions are not so concentrated 
near the origin. 


The effect of the tensor interaction Eq. (92) is very small in the case of the hC 


wave functions, so there is no way to determine the parameter B directly from 
the spectrum. However, assuming for B a value of about 1/10 A, the description 
of the spectrum is slightly improved. The final results for the spectrum are 
reported in Fig. ([^. 

The description of the spectrum is particularly good for the negative parity 
resonances. The analytical model has been also used for the calculation, without 
free parameters, of the electromagnetic transition amplitudes to some negative 
parity baryon resonances [531 HI] ■ 

The success of the 1/x model in describing the spectrum indicates that, at 
least for the inner states, the confinement is provided by the hC potential [94] . 
as further supported by ref. [55] . 


4 The baryon spectrum 

4.1 Results from the hCQM 

As discussed in the previous section, the free parameters of the hCQM are 
determined by fitting the masses of the 4* and 3* resonances m reported in 
Fig. (0. Of course the model can predict the masses of all other resonances 
belonging to the first three energy shells and the number of theoretical states 
exceeds the observed one, leading to the problem of the missing resonances, 
a problem in common with other CQM, in particular the h.o. one. In this 
respect, it is interesting to compare the number of states predicted by the two 
potentials, the h.o. (see Table0 and the hCQM, which, as shown in Table|^are 
30 and 39 respectively. These two values are certainly larger than the number 
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of observed 4* and 3* states, however the situation becomes different if one 
considers separately the positive and negative parity states and if also the new 
results from PDG2012 [63] are taken into account. 


Table 3: The number of states predicted by the h.o. and hCQM models, re¬ 
porting separately the positive and negative parity N and A states. In the last 
four columns, the number of states listed in the 2010 |6T| and 2012 [^ PDG 
editions are reported. 



h.o. 

hCQM 

PDGIO 
4* -t 3* 

PDGIO 

4*-t3*-k2* 

PDG12 
4* -t 3* 

PDG12 

4*-t3*-k2* 

N+ 

14 

15 

5 

8 

6 

11 

N- 

5 

10 

5 

7 

6 

9 

A+ 

9 

10 

6 

7 

6 

7 

A- 

2 

4 

2 

3 

2 

4 

Total 

30 

39 

18 

25 

20 

31 


The positive parity states allowed by the h.o. model are abundant, but the 
negative ones are just what is necessary for the description of the observed 4* 
and 3* states of PDG2010 [ST]. In the hCQM there are two positive parity 
states more and the negative parity ones are doubled, in agreement with the 
fact that in the hGQM spectrum (see Fig. 14) there are two extra levels, one 0^ 
with a Pll and a P33 state, and a level, where a further series of negative 
parity states can be settled. 

In the last edition of the PDG [S3], some new states are reported. This 
achievement has been possible also thanks to the availability of very precise cross 
section and polarization data from photoproduction experiments at GLAS m 
and by the most recent coupled- channel analysis of the Bonn-Gatchina group 
|101| (for a discussion see |102p . In particular there is a resonance IIi3(1875) 
with a 3* status. The negative parity states allowed by the h.o. model are all 
already occupied and only the hGQM, with its further negative parity level, can 
describe such a state. Moreover, if one considers also the new 2* states [H3] . 
there are 13 negative parity resonances, 9 of the N and 4 of the A type, to be 
compared with the allowed values of the hCQM, that is 10 and 4, respectively. 
Finally, the total number of observed 4*, 3* and 2* resonances is 31, greater 
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than the number allowed by the h.o. and not so far from 39, the hCQM value. 

The comparison of the theoretical spectrum with all the 4*, 3* and 2* listed 
in PDG [63] is shown in Fig. ([^ |103j . In this Figure two theoretical levels are 
not shown, that is one A^l/2 and one N3/2 state belonging to the (20,1+) SU (6) 
multiplet; they are mixed by the hyperfine interaction only with the states with 
7 + z/ greater than 2 and do not contribute to the structure of the nucleon. 


(j0V *** 



Figure 9: (Color online) The spectrum obtained with the hCQM described in 
Sec. 3 in comparison with all the 4*, 3* and 2* resonances reported in the 2012 
edition of the PDG [55] . 

The overall description of the spectrum is quite good, considering that the 
model has only three free parameters. As for the Roper resonance, it is practi¬ 
cally degenerate with the negative parity states, thanks to the behaviour of the 
hC potential, which is only slighlty modified by the confinement term. However, 
the theoretical Roper mass is still to high with respect to the experimental data. 
This problem is common to various CQM and it has some time ago suggested 
the idea that the Roper is not simply a “breathing” mode of the nucleon, but 
it is a hybrid state qqqG [104] . that is three quark plus a gluon component. 
However this model for the Roper has been ruled out by the recent results on 
the 7 —7> iV(1440) data |105j . which showed that the longitudinal electroex¬ 
citation is significantly non zero while the hybrid model predicts that it should 
be vanishing. 

The theoretical states in the higher part of the spectrum are somewhat 
compressed as an effect of the hC interaction. However, the masses of the 
resonances in this region are determined with large uncertainties and, because 
of their strong decays, the states are expected to have large widths and to be 
partially overlapped. When comparing the theoretical baryon spectra with the 
experimental data one should not forget that up to now GQM models predict 
states with zero width, since no coupling to the continuum has been consistently 
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introduced. Some time ago the Isgur-Karl model has been implemented with 
quark-meson couplings [THH] . allowing to calculate both the strong decay widths 
and the effects of the continuum on the resonance energies. This approach 
considers quark and meson degrees of freedom on the same footing, it would be 
desirable on the contrary to have an approach in which quark-antiquark pair 
mechanisms are consistently taken into account. In this respect an important 
improvement has been achieved by a recent work [ina [Ml [Ml, in which an 
unquenched constituent quark model for baryons has been formulated and the 
quark-antiquark pair contributions are taken into account consistently. 


4.2 The hCQM with isospin 

In the hCQM hamiltonian, the SU{6) violating term is provided by the hyper- 


fine interaction of Eq. (19), similarly to what happens with the Isgur-Karl model 

and its semirelativistic extension |7]. However the SU{6) violation can be 
given also by a flavour dependent term, which is more or less explicitly included 
in other approaches. In fact, in the algebraic model (BIL [8]) the quark energy 
is written in terms of Casimir operators of symmetry groups which are relevant 
for the three-quark dynamics (see Sec. 2.4). In particular, for the internal de¬ 
grees of freedom the Giirsey-Radicati mass formula m is used, leading to an 
isospin dependent term which turns out to be important for the description of 
the spectrum. In the yCQM [12], the quark-quark interaction is provided by 
one meson exchange and therefore the corresponding potential is spin-flavour 
dependent and is crucial for the description of baryons up to 1.7 GeV. As for 
the BN [T7]) model, the SU{6) violation arises from the instanton interaction 
which does not act on A states. Moreover, it has been pointed out that an 
isospin dependence of the quark potential can be obtained by means of quark 
exchange m- 

Therefore there are many motivations for the introduction of a flavour de¬ 
pendent term in the three-quark interaction and for this reason also in the case 
of the hCQM an isospin dependent term has been included in the quark inter¬ 
action [ml dm. 

The SU(6) violation coming from the hyperfine interaction is still present, 
with one important modification, namely in the spin-spin interaction the 5-like 
term is substituted with a smearing factor given by a gaussian function of the 
quark pair relative distance m- 


Hs 


As E 

i<j 


{y/nasf ® 


(s' -Sj). 


(93) 


The remaining SU{6) violation comes from two terms. The first one is 
isospin dependent 


Hi 


i<j 


1 

iV^air 



(94) 
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where ti is the isospin operator of the i-th quark and rij is the relative quark 
pair coordinate. The second one is a spin-isospin interaction, given by 


Hsi = Asi 


i<j 




e "SI (si ■ Sj){ti -tj), 


(95) 


where s) and U are respectively the spin and isospin operators of the i-th quark 
and rij is the relative quark pair coordinate. The complete interaction is then 
given by 

H,nt = Vix)+Hs + Hi + Hsi, (96) 


where V(x) is the hypercentral potential of Eq. (811, 
for the 3*- and 4*- resonances is shown in Fig. (10). 


The resulting spectrum 



Figure 10: (Color online) The spectrum obtained with the hCQM with the 
spin and isospin dependent interactions of Eq. (961 in comparison with all the 
4*, 3* and 2* non strange resonances reported in the 2012 edition of the PDG 


The TV — A mass difference is no more due only to the hyperfine interaction. 
In fact, in this model its contribution is only about 35%, the remaining splitting 
coming from the spin-isospin term (50%) and from the isospin one (15%). 

It should also be noted that the negative parity resonances are again well 
described. In this model however there is the correct inversion between the 
Roper and the negative parity resonances and this is almost entirely due to the 
spin-isospin interaction, as stated in Ref. m- In general, the position of the 
Roper resonance is reproduced in all models containing an isospin dependent 
interaction ismiiiT]. 

Also the higher states are fairly described and slightly less compressed than 
in the standard hCQM. 
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The tensor part of the hyperfine interaction, which is omitted for simplicity 
in Eq. (96), is taken into account in the calculation, however its contribution to 
the spectrum is negligible. 


4.3 An extension to strange baryons 


The hypercentral interaction of Eq. ( |81[ ) describes the average energies of the 
SU{6) multiplets, while the splittings within each multiplets are generated by 
the hyperfine interaction Eq. (19) or by the spin-isospin interaction of Eq. (96). 
In the latter case the flavour dependence is due only to the isospin operators, 
provided that the interest is limited to the non strange baryons. In order to 
describe the spectrum of strange baryons as well, it is necessary to introduce a 
flavour dependence which involves both isospin and strangeness. This can be 
achieved in the hCQM in a similar manner to the algebraic model [8] quoted in 
Sec. 2.4, that is describing the SU{6) violation by means of a Giirsey-Radicati 
(GR) mass formula |I13j . 

The original GR mass formula m can be rewritten in terms of Gasimir 
operators m 


M = Mo+ C C2[SUsi2)] + D C'i[C/y(l)] -h E 


C2[SUi{2)] 


i(Ci[C/y(l)])2 


(97) 

where C2[SUs{2)] and C2[SUi{2)] are the SU{2) (quadratic) Gasimir operators 
for spin and isospin, respectively, C'i[?7y(I)] is the Gasimir for the U{1) subgroup 
generated by the hypercharge Y. 

However, in the framework of the GQM, the underlying symmetry is provided 
by SU{6) and Eq. ( [9^ is not the most general formula that can be written on 
the basis of a broken SU{6) symmetry. It can then be generalized as follows 

m 


M = Mo +4 C2[SUsf{6)] + B C2[^Gj’(3)] + C C2[SUs{2)] + 

+D Ci[C/y(l)] + E {C2[SUi{2)] - i(Ci[t/y(l)])2), (98) 

where Mq is the SUspi^i) invariant mass. 

The idea is then to consider the energy levels provided by the hypercentral 
potential of Eq. ( [M] ) as the values of the central masses of the SUsf{G) mul¬ 
tiplets and to use the generalized mass formula in order to describe the spin- 
flavour splittings within the multiplets m- The hamiltonian is assumed to 
be 

H = Ho + Hgr, (99) 

where iJo is the hCQM hamiltonian without the hyperfine interaction 

—* 2 2 

Elo = 3m + -I- ^ - — + ax (100) 

2m 2m x 

and ELgr is given by the spin-flavour dependent part of Eq. ( [M| ) 

Hgr = a C 2 [SUrf (6)] + B C2 [^C/f (3)] + C C2 [SUg (2)] + 

YD Ci[Uy{1)]+E {C2[SUi{2)] -\{Ci[Uy{IW) ■ (101) 
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Table 4: The eigenvalues of the quadratic Casimir operators for the groups 
SlIsriG) (left) and SUf{^) (right). 


SUsf{Q) 

C 2 

S'Pf(3) 

C 2 

[56] 

45 

4 

[8] 

3 

[70] 

33 

4 

[10] 

6 

[20] 

21 

4 

[1] 

0 


In order to apply the generalized GR mass formula to the baryon spectrum 


it is necessary to assume that the coefficients A,B, .. .in Eq. (101) be the same 


in the various SUsf{(^) multiplets. This actually seems to be the case, as shown 
by the algebraic approach to the baryon spectrum [8], where a formula similar 
to Eq. (llOll) has been applied. 


The matrix elements of Hqf are completely determined by the values of the 
various Casimir operators |113j : for the SUsf{Q) and SUf{5) groups the values 
of the Casimir operator C 2 are reported in Table while for the SU{2) and 
C/y(l) groups one has 

{C2[SUi(2)])=I{I+1), (Ci[[/y(l)])=y, (C2[5 C/s(2)]) = ,5(,5+1). (102) 


The mass of each baryon state is then 


{B\H\B)=E^, + {B\Hgr\B), 


(103) 


where E^i, are the eigenvalues of the hypercentral potential of Eq. (100). 


The parameters a and t of the hypercentral potential have been fitted in 
Sec. 3.3 in presence of the hyperfine interaction. Here the SU{6) violation is 
provided by a different mechanism and then these parameters must fitted to the 


spectrum together with those introduced in Eq. (101). 


Such fit can be performed in two ways. The hrst one is an analytical proce¬ 
dure which consists in choosing a limited number of well known resonances and 
expressing their mass differences in terms of the Casimir operator values. In this 
way a part of the unknown coefficients is evaluated directly, while the remaining 
ones is fitted to the experimental spectrum. A possible choice of resonance pairs 
is given by the following ones 


(fV(1650)S'll - 
(A(1232)P11 
(A(1535)S'll - 
(E(1193)P11 
(A(1116)P01 


- A(1535)S'll) 

- 7V(938)P11) 
7V(1440)P11) 

- 7V(938)P11) 

- iV(938)Pll) 


3C, 

9B + 3C + 3P, 
Eio — Eqi + 12A, 
3/2P- P, 
-D-1/2E. 


(104) 
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Figure 11: (Color online) The spectrum obtained with the hCQM and the 
generalized Giirsey-Radicati formula of Eq. (981 in comparison with all the 4* 
and 3* N, A and A resonances reported in the 2000 edition of the PDG |114] . 
The parameters given by the direct fit are used (see Table column (11)). 


I I 



Figure 12: 


(Color online) The same as in Fig. 11 for the E, 2 ; and resonances. 
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Table 5: The fitted values of the parameters of the Hamiltonian (99) |113j . 
Column (I) reports the values given by the analytical procedure of Eq. (104), 
while Column (II) contains the results of the direct fit to the experimental 
spectrum 


Parameter 

(I) 

(11) 

a 

I.4/m“^ 

2 .1 /to“^ 

T 

4.8 

3.9 

A 

-13.8 MeV 

-11.9 MeV 

B 

7.1 MeV MeV 

11.7 MeV 

C 

38.3 MeV 

30.8 MeV 

D 

-197.3 MeV 

-197.3 MeV 

E 

38.5 MeV 

38.5 MeV 


In the second procedure all the parameters are fitted in order to reproduce 
the baryon spectrum. The resulting values of the parameters are reported in 
Table and the corresponding spectrum is shown in Fig. El The overall de¬ 
scription of the baryon spectrum of the 4* and 3* resonances is quite good, 
specially considering the simplicity of the model. 

Using the analytical procedure the overall agreement with the spectrum is 
slightly worsened, but the strange sector is better described. 

In both procedure, there is the need of a non zero value of the parameter 
A in order to reproduce the spectrum. An attempt to fit the data fixing A=0 
has been tried, however the resulting parameters a and t are quite different 
form those reported in Table Furthermore, the correct ordering of the Roper 
resonance and the negative parity resonances is lost. The presence of the Casimir 
C' 2 [S'Us_f( 6 )] is essential in order to shift down the energy of the first excited 
O’*" state with respect to the I“, an effect which is similar to the one produced 
by the U-potential in the Isgur-Karl model and by the presence of a flavour 
dependent interaction in the previously quoted CQMs. 


5 The electromagnetic excitation of baryon res¬ 
onances 

5.1 The transition amplitudes 

The electromagnetic excitation of the baryon resonances is an important source 
of information concerning the nucleon structure. The absorption of real photons 
is a direct measure of the excitation strength while the inelastic electron scat- 
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tering is a probe of the excited nucleon structure at short distances. There are 
presently many experimental data taken at various laboratories (Jlab, Mainz, 
Bonn, ...) but a systematic study of the electromagnetic excitation of the nu¬ 
cleon at high is expected to be performed by the upgraded 12 GEV beam 
at Jlab mnuni]. 

There is an intense theoretical and phenomenological activity which aims 
at extracting the transition amplitudes from the experimental data on photo- 
and electro-production of mesons off nucleons using mainly the Partial Wave 
Analysis m- Various groups have devoted much effort in this sense, using 
different techniques. Among them we quote the analyses made by the following 
groups: George Washington University (SAID) |118j . Mainz University (MAID) 
|119LII20l 1121) . Dubna-Mainz-Taipei (DMT) |I22l 1120] Bonn-Gatchina (BnGa) 
EBAG at Jefferson Lab Jiilich Giessen Zagreb-Tuzla 

m- 

Erom the theoretical point of view, the photo- and electro-excitations of the 
nucleon to the various baryon resonances are described by the helicity ampli¬ 
tudes, defined as the matrix elements of the electromagnetic interaction, 
between the nucleon, N, and the resonance, B, states: 

Ai/ 2 = y^(i3,J',J'= i |J+|iV,J = = -i ), 

A 3/2 = = i |J+|iV,J = i,J, = i), (105) 

5i/ 2 = ^(i3,J',J' = i |Jo|iV,J = i,J. = i), 


where ko is the photon energy and, for the transverse excitation, the photon has 
been assumed, without loss of generality, as left-handed. 

The hGQM has been completely specified in Sec. 3.3. We recall that the 


Hamiltonian is given by Eq. (82) 


HhCQM = Sm+ ^ -I- ^ 
2 m 2m 


- \- ax - 


H, 


hyp, 


(106) 


where a = 1.16 fm“^, r = 4.59 and the strength of the hyperfine interaction is 
fixed by the N — A mass difference. 

The states of the various resonances have been explicitly built up and there¬ 
fore they can be used for the calculation of any quantity of physical interest. In 
order to proceed to the calculation of the helicity amplitudes, one has to specify 
the current in the electromagnetic interaction. In the framework of the hCQM, 
the current is simply given by the sum of the quark currents j^{i) 




( 107 ) 
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and will be used in its nonrelativistic form 


p{k) = 


2 = 1 


1 


3 


2 = 1 


ia^{i) X k)e^ 


where is the charge of the i-th quark 

2^3 


Si = 


+ 'r|(*)]: 


(108) 


(109) 


f) is the quark coordinate and are, respectively, the quark spin and 

isospin operators. 

In order to compare the theoretical results with the experimental data, 
the calculation should be performed in the rest frame of the resonance (see 
e.g. [HH]). The nucleon and resonance wave functions are actually calculated 
in their respective rest frames and, before evaluating the matrix elements given 


in Eqs. (105), one should boost the nucleon to the resonance c.m.s.. However, 


in order to minimize the discrepancy between the nonrelativistic and the rel¬ 
ativistic boosts in comparing with the experimental data, we can consider the 
Breit frame, as in refs. ( [33 |H1 SS] ). In this frame = —PR = —fc/2, where 
Pn, Pr and k are, respectively, the nucleon, resonance and photon trimomenta. 
The relation of the latter with the momentum transfer squared is given by: 


k^ = Q^ + 


2(M2 -h IT2) -h Q2 


( 110 ) 


— fi,g. 


where M is the nucleon mass, W is the mass of the resonance and 
ko being the photon energy. 

Furthermore, one has to consider that the helicity amplitudes extracted from 


the photoproduction contain also the sign of the ttNN* vertex (see Fig. 13). The 
theoretical helicity amplitudes are then defined up to a common phase factor ( 


^1/2,3/2 — C -^1/2,3/2 ^ 1/2 — C ^ 1 ^ 2 - (HI) 

The factor ^ can be taken |3S] in agreement with the choice of ref. [32] , with the 
exception of the Roper resonance, in which case the sign is in agreement with 
the analysis performed in [129] . 

In the following Subsections the results for the photo- and electro-excitation 
of the baryon resonances calculated with the hCQM will be presented. For 
consistency reasons, in the calculations the values of W given by the model are 
used instead of the phenomenological ones. 

The calculations regarding the photocouplings have already been published 
in [23], those for the transverse excitation of the negative parity resonances in 
[33] j while a systematic study of all the electromagnetic excitations has been 
reported in [3S|- 

It should be stressed that the results are obtained with a parameter free 
calculation, that is they are predictions of the model. 
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Table 6: Photocouplings (in units 10“^ GeV“^/^) predicted by the hCMQ [46] 
in comparison with PDG data [63] for proton excitation to N*-like resonances. 
The proton transitions to the 811(1650), D15(1675) and D13(1700) resonances 
vanish in the SU{6) limit. 


Resonance 

"^1/2 

hGQM 

^1/2 

PDG 

AP 

^3/2 

hCQM 

A'P 

PDG 

Pll(1440) 

88 

-65 ±4 



D13(1520) 

-66 

-24 ±9 

67 

166 ±5 

811(1535) 

109 

90 ±30 



811(1650) 

69 

53± 16 



D15(1675) 

1 

19 ±8 

2 

15 ±9 

F15(1680) 

-35 

-15 ±6 

24 

133 ± 12 

D13(1700) 

8 

-18± 13 

-11 

-2 ±24 

Pll(1710) 

43 

9 ±22 



P13(1720) 

94 

18 ±30 

-17 

-19 ±20 


Table 7: The same as in Table [^ but for the N — A excitation. 


Resonance 

"^1/2 

hCQM 

PDG 

AP 

^3/2 

hCQM 

AP 

^3/2 

PDG 

P33(1232) 

-97 

-135 ±6 

-169 

-250 ±8 

831(1620) 

30 

27± 11 



D33(1700) 

81 

104 ±5 

70 

85 ±2 

F35(1905) 

-17 

26± 11 

-51 

-45 ± 20 

F37(1950) 

-28 

-76 ± 12 

-35 

-97± 10 
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^ 1/2 ^ 3/2 ^ 1/2 


Figure 13: (Color online) The electromagnetic excitation of nucleon resonances. 
The decay can occur in more than one pion. 


5.2 The photocouplings in the hCQM 

The resonances which have been considered are those which, according to the 
PDG classification [63], have an electromagnetic decay with a three- or four- star 
status. This happens for twelve resonances, namely the positive parity I = \ 
states 

Pll(1440), F15(1680), Pll(1710), (112) 

the negative parity I = \ states 

1113(1520), S'll(1535), S'lhlGSO), 1115(1675) (113) 

and the 7 = | ones 

P33(1232), S'31(1620), 1133(1700), F35(1905), F37(1950). (114) 


Besides these states, we have considered also the two resonances D13(1700) and 
P13(1720), which are excited in an energy range particularly interesting for the 
phenomenological analysis. 

The proton and neutron photocouplings predicted by the hCQM |24| are 
reported in Tables EE and in comparison with the PDG data [S3] and the 
analysis by the Bonn-Gatchina group m- The overall behaviour is fairly well 
reproduced, but in general there is a lack of strength. The proton transitions to 
the 811(1650), D15(1675) and D13(1700) resonances vanish exactly in absence 
of hyperfine mixing and are therefore entirely due to the SU (6) violation. 

As already noted in Sec. 2.2, the hyperfine interaction is responsible for a 
deformation of the A resonance and therefore the ratio of Eq. (231 is different 
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Table 8: Photocouplings (in units 10“^ GeV“^/^) predicted by the hCMQ [3^ 
in comparison with PDG data [63] and the recent Bonn-Gatchina analysis |130j 
for neutron excitation to N*-like resonances. 


Resonance 

An 

^1/2 

hCQM 

An 

^1/2 

PDG 

An 

^1/2 

BnGa 

An 

^ 3/2 

hCQM 

An 

^ 3/2 

PDG 

An 

^ 3/2 

BnGa 

PI1(1440) 

58 

40± 10 

43± 12 




D13(1520) 

-1 

-59 ±9 

-49 ±8 

-61 

-139 ± 11 

-113 ± 12 

811(1535) 

-82 

-46 ± 27 

-93± 11 




811(1650) 

-21 

-15 ±21 

-25 ±20 




D15(1675) 

-37 

-43 ± 12 

-60 ±7 

-51 

-58± 13 

-88± 10 

F15(1680) 

38 

29± 10 

34 ±6 

15 

-33 ±9 

-44 ±9 

D13(1700) 

12 

0±50 


70 

-3 ±44 


Pll(1710) 

-22 

-2± 14 

-40 ± 20 




P13(1720) 

-48 

1±15 

-80 ±50 

4 

-29 ±61 

-44 ±9 


from zero. This ratio can also be expressed in terms of the helicity amplitudes 


Rem 


Ge _ 's/s Ai/2 — ^3/2 
Gm Vs ^ 1/2 + 3 ^ 3/2 


(115) 


with the theoretical values reported in Tablej^ Rem turns out to be smaller than 
the experimental one. The point is that the E2 transition strength predicted by 
the hCQM is too low and a possible explanation of this result will be discussed 
later. 

The results obtained with other calculations are qualitatively not much dif¬ 
ferent (231 US] and this is because the various GQM models have the same SU (6) 
structure in common. 

It should be reminded that in previous nonrelativistic calculations with h.o. 
wave functions m, it was necessary to assume a proton radius of the order of 
0.5 fm in order to ensure a vanishing ^^^2 resonances D13(1520) and 

F15(1680), whose peaks are absent in the forward photoproduction [21]. The 
proton radius calculated with the hGQM is actually 0.48 fm and this explains 
why the predictions of the hGQM do not differ too much from the other calcu¬ 
lations. A too low proton radius is of course a problem if one wants to calculate 
the elastic form factors of the nucleon, but for the description of the helicity 
amplitudes it is beneficial and, as we shall see later, it plays an important role 
in the discussion concerning the mechanisms which are missing in any CQM. 

There are now many new analyses concerning the neutron helicity amplitudes 


41 









(' |130j and references quoted therein). In Table we report also the results of 
the Bn-Ga analysis USD]. The hCQM predictions are in fair agreement with 
these data, perhaps better than with the PDG ones. 



(GeV^) 


▲ 

(c) hCQM - 

- 

Maid07 ^ 

. 

Azn09 

4 






0 1 2 3 4 5 

(GeV^) 



Figure 14: (Golor on line) The P33(1232) helicity amplitudes predicted by the 
hCQM (full curves) A 3/2 (a), ^ 1/2 (b) and S' 1/2 (c), in comparison with the data 
of ref. |131] and with the Maid2007 analysis |121j of the data by refs. |132j and 
m- The PDG points [53] are also shown. The figure is taken from ref. [35] 
(Gopyright (2012) by the American Physical Society) 


5.3 The helicity amplitudes in the hCQM 

As already mentioned, a systematic review of the hGQM predictions for the 
transverse and longitudinal helicity amplitudes and their comparison with the 
experimental data is reported in ref. |4tij . Here we limit ourselves to some of 
the most important excitations. 

The transition amplitudes for the excitation of the P33(1232) resonance are 
given in Fig. 

The transverse excitation to the A resonance has a lack of strength at low Q^, 
a feature in common with all GQM calculations. The medium-high behavior 
is decreasing too slowly with respect to data, similarly to what happens for the 
nucleon elastic form factors [301132]. As we shall see later, the nonrelativistic 
calculations are improved by taking into account relativistic effects. Since the 
A resonance and the nucleon are in the ground state S'[/(6)-configuration, we 


42 




















(GeV^) (GeV^) 


Figure 15: (Color on line) The Pll(1440) proton transverse (a) and longitndinal 
(b) helicity amplitudes predicted by the hCQM (full curves), in comparison with 
the data of refs. m, [HI] and the Maid2007 analysis |121) of the data by 
refs. |135j . |132] . |133) and |136j . The PDG point (53] is also shown. The figure 
is taken from ref. [23] (Copyright (2012) by the American Physical Society). 


expect that their internal structures have strong similarities and that a good 
description of the N—A transition from factors is possible only with a relativistic 
approach. Such feature is further supported by the fact that the transitions to 
the higher resonances are only slightly affected by relativistic effects [30] . 

The Roper excitation is reported in Fig. 15 Because of the - term in the 


hypercentral potential of Eq. (82), the Roper resonance can be included in the 
first resonance region, at variance with h.o. models, which predict it to be a 2 hui 
state. There are problems in the low region, but for the rest the agreement 
is interesting, specially if one remembers that the curves are predictions and the 
Roper has been often been considered a crucial state, non easily included into a 
constituent quark model description. In particular, the longitudinal excitation 
is quite different from zero nna, in agreement with the hCQM and at variance 
with the hybrid qqq-gluon model [104j . In the present model, the Roper is a 
hyperradial excitation of the nucleon. 

We consider now the excitations to some negative resonances |2SJ|25|, namely 
the D13(1520) and the 811(1525) ones, reported in Figs. 16 and 17 respectively. 

The agreement in the case of the Sll is remarkable, specially if one considers 
that the hCQM curve for the transverse transition has been published three 
years in advance |45| with respect to the recent TJNAF data |131j . |139j . |141] . 


It is interesting to discuss the influence of the hyperfine mixing on the ex¬ 
citation of the resonances. Usually there is only a small difference between the 
values calculated with or without hyperfine interaction. In some cases, however 
the excitation strength vanishes in the SU(6) limit, as already mentioned in 
Table [^ the non vanishing result is then entirely due to the hyperfine mixing 
of states. In the case of the 811(1650) resonance, the resulting transverse and 
longitudinal excitations have a relevant strength. 

The three helicity amplitudes of the D13(1700) resonance are again non zero 
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Figure 16: (Color on line) The D13(1520) proton helicity amplitudes predicted 
by the hCQM (full curves) A 3/2 (a), Aij^ (b) and S' 1/2 (c), in comparison with 
the data of refs. [Hi], [m], with the compilation reported in refs. [13711138] 
and the Maid2007 analysis |121| of the data by refs. [132] and [133] . The PDG 
points l^are also shown. The figure is taken from ref. (JS) (Copyright (2012) 
by the American Physical Society) 


because of the hyperfine mixing but the excitation strength is very low. Also 
in the case of the transverse excitation of the D15(1675), the strength is given 
by the hyperfine mixing, while the longitudinal amplitude iS' 1/2 vanishes also in 
presence of a SU ( 6 ) violation. 


The results of the hCQM, as reported in Figs. 14 15 16 17 and in ref. 


show some general features. First of all, there is a nice agreement with data at 
medium-high values of . This is mainly due to the presence of the l/cc term 
in the hCQM hamiltonian. In fact, the helicity amplitudes calculated with the 
analytical model of refs. |S31|Mj have a dependence very similar to the com¬ 
plete hCQM and numerical values only slightly different. The wave functions of 
the analytical model are exactly those determined by the 1 /x potential, which 
gives then the main contribution to the transition strength. 

Another important fact which helps in obtaining a good behaviour of the 
hCQM helicity amplitudes is the smallness of the resulting proton radius. As 
already mentioned, the r.m.s. radius of the proton calculated with the hCQM 
wave functions corresponding to the parameters of Eq. (85) is 0.48 fm, which is 
very near to the value necessary in order to fit the D13 photocoupling [STj . 
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Figure 17: (Color on line) The 811(1525) proton transverse (a) and longitudinal 
(b) helicity amplitudes predicted by the hCQM (full curve), in comparison with 
the data of refs. |139j (open diamonds), |131j (full diamonds), |140j (crosses), 
|141j (open squares), |142) (full squares), the Maid2007 analysis | 121 ) (full tri¬ 
angles) of the data by refs. |132j and the compilation of the Bonn-Mainz-DESY 
data of refs. [14311144(114511146] (stars), presented in |142j . The PDG point [HH] 
(pentagon) is also shown. The figure is taken from ref. [46] (Copyright (2012) 
by the American Physical Society) 


Both these features, the presence of the hyperCoulomb term in the quark 
potential and the smallness of the proton radius, concur in order to obtain the 
results shown in the figures. This can be seen also looking at Fig. where 
the results of Fig. ([T^ for the D13(1520) excitation are given in comparison 
with the h.o. curves corresponding to two values of the proton radius, namely 
the experimental one (about 0.86 fm) and the one fitted to the amplitude 
(about 0.5 fm). The curves with the correct proton radius are completely out 

of the experimental data because of the gaussian factor e soS ^ typical of the 
h.o.; a is the h.o. constant (see Sec. 2.2) and its value corresponding to the 
experimental proton radius is 0.229 GeV. In the case of the smaller radius 0.5 
fm, with a = 0.41 GeV, the ^ 3^2 amplitude is of course well reproduced but the 
h.o. curve for the A^y 2 amplitude is again far from the data. For this reason, 
we expect that the hCQM should be a good starting point for the description 
of the non perturbative components of the parton distribution. 

There are in general discrepancies at low displaying a lack of strength 
which is typical of all GQMs. Nevertheless, in many cases the ^^^2 amplitudes 
are better reproduced than the ^ 3^2 ones. 

These shortcomings of the hCQM results could be ascribed to the non- 
relativistic character of the model. In fact, the electromagnetic excitation leads 
to a recoil of both the nucleon and the resonance and the effect is expected to 
increase with the momentum transfer while the wave functions are calcu¬ 
lated in the rest frame of each three-quark system. As it will be discussed in 
Sec. 6.2, it is possible to apply Lorentz boosts in order to bring the nucleon and 
the resonance to a common Breit frame, but this relativistic corrections produce 
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Figure 18: (Color on line). The hCQM predictions for the transverse D13(1520) 
helicity amplitudes (full curve) in comparison with the h.o. results corresponding 
to two values of the proton radius: 0.5 fm (hoi) and 0.86 (ho2). The data are 
the same as in Fig. 


only a slight modification of the hCQM results 

There is a consensus on the fact that the missing strength at low is due to 
the lack of quark-antiquark effects [45] , probably important in the outer region 
of the nucleon. This statement receives a strong support by the explicit calcula¬ 
tions of the meson cloud contributions to the helicity amplitudes performed in 
the framework of dynamical models (see ref. [12011148] and references therein). 

In particular the Dubna-Taipei-Mainz (DMT) model introduces the pion 
cloud contribution to the electromagnetic excitation according to the mechanism 
shown in the upper left part of Fig. (191. In the same hgure, the longitudinal 
Si /2 and transverse ^ 3/2 and A 3/2 helicity amplitudes for the TV—A are reported 
m- The theoretical predictions of the hCQM (full curves) are compared with 
the results of the MAID fit and with the pion cloud contribution calculated 
by means of the DMT model (dashed curve) [120] . The hCQM results are 
much lower than the experimental data, but the pion cloud contribution gives 
relevant contributions just where the hCQM is lacking. This is particularly 
evident in the case of the longitudinal iS' 1/2 amplitude: the hCQM predicts an 
almost vanishing value while the pion alone seems to be able to account for 
the data. Of course one cannot simply add the hCQM and pion contributions, 
since they are calculated in two different and inconsistent frameworks, but it is 
nevertheless interesting that the pion cloud seems to contribute systematically 
where the hCQM is lacking, as it can be seen also for the helicity amplitudes of 
many other resonances ma¬ 
in this way the emerging picture in connection with the electromagnetic 
excitation of the nucleon resonances is that of a small confinement zone of about 
0.5 fm surrounded by a quark-antiquark (or meson) cloud. The calculations 
of the meson cloud performed with the DMT model shows that this picture 
seems to be reasonable, but the problem is how to include in a consistent way 
the quark-antiquark pair creation mechanisms in the framework of the CQM. 
This goal can be achieved by unquenching the quark model and, as already 


46 














-pion cloud contribution 






(GeV/c)’ 



Q» (GeV/c)’ 



Figure 19: (Color on line). The pion cloud mechanism considered in the DMT 
model |120j (upper left). The remaining figures report the hCQM predictions 
for the TV — A excitation (full curves) in comparison with the meson cloud 
contribution (dashed curves) and the results of the MAID fit |150| (full points 
and the htted curves passing through them) |149j . 


mentioned earlier, an important improvement has been achieved by a recent 
work [1071110811109j . We shall come back on this point in the discussion. 

6 The elastic form factors of the nucleon 

6.1 Introductory remarks 

An important aspect of the quark model predictions concerning the elastic form 
factors is their dependence, which is strictly related to the form of the quark 
wave functions and then of the quark potential. 

We can start studying the nucleon charge form factor in absence of the 
hyperfine mixing. The nucleon state (see App. B) can be written as 

\N) = '0oo(a:)D[o]oo^(XMS<(>MS + Xmai/^ma), (116) 

where D[o]oo = l/( 47 r) 24 /v^. 
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The nucleon charge form factor is given by the matrix element of the charge 


density operator of Eq. (108) 


= {N\3 


1 + 3To(3) 


6 


1^), 


(117) 


where to( 3) is the third component of the isospin operator of the third quark, 
and in the Breit system, according to Eq. (llOl. Introducing 


the hyperspherical coordinates and performing the integrals over the angle and 
hyperangle variables, one gets 




1 + To 


F(fc), 


where tq is the third component of the nucleon isospin and 


m = 




dx ijjooix)'^ J 2 {o,kx), 


(118) 


(119) 


J 2 {z) being a Bessel function of integer order 2. Eq. (119) can be inverted 


obtaining an expression of the wave function in terms of the form factor 


^4 rao 

tpooix)'^ = ^ J dk k^ F{k) J2{akx). 


( 120 ) 


In this way, starting from any given form factor it is possible to obtain the 
appropriate wave function iPqq{xY and then also the potential for which ipooix) 
is the ground state. Assuming the dipole form factor 1/(1 + b'^k'^), the resulting 
potential is given by 


a",. l.A:o(y).2 , Ko{y) 


^‘'*^(^) “ 2&2 [1 2^K^{y) 


Y-A 


6 


yKi{y) y"^ 


( 121 ) 


where y = ax/h and Ko{y), Ki{y) are the modified Bessel function of the second 
kind. It is interesting to note that for large values of y, the potential assumes 
the form 


Vd. 


tp 


“ n ^ 1 

(^) Tl9 - 9 ] 

41)2 y y2 J 


( 122 ) 


that is there is no confinement. 

The nucleon charge form factor turns out to be proportional to the charge 
and therefore it is zero for the neutron. This happens as long as the space part 
of the state is completely symmetric. Because of the hyperfine interaction, also 
the state state \N'^Sif 2 )m, having mixed space symmetry, contributes to the 
nucleon (see Eq. (83)), thereby generating a non-zero neutron form factor |68j . 


In h.o. models the ground state wave function is given by a gaussian e 2 

_ 

which leads to the form factor e ei^, where a is the h.o. constant (see Sec. 2.2). 
Because of the hyperfine mixing, the nucleon state is given by a superposition 


of SU (6) configurations of Eq. (83), but the dominant behaviour of the form 
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factor is in any case a gaussian, which is too strongly damped with respect to 
the experimental data. 

In the case of the purely hypercoulomb potential, the ground state wave 
function is given by an exponential function and the corresponding form 
factor is 

Gl(.Q^)\hC = I 2 , ■ (123) 

The use of the hCQM with the confinement potential modifies this power-law 
behaviour, however, even taking into account the hyperfine mixing, the resulting 
form factor has a more realistic behaviour with respect to the h.o. model. 

As for the other nucleon form factors, in the SU{6) limit, there is perfect 
scaling, similarly to the dipole fit, in the sense that 

G'^(Q^)|/iCQm = —Gm(*3^)Uc’qm = — G^{Q‘^)\hCQM (124) 

f^p l-^n 


and 

Gl{Q^)\hCQM = 0. (125) 

The hyperfine mixing modihes scarcely the above results, with the already 
mentioned difference of producing a non zero neutron charge form factor. 

Comparing the predictions of the hCQM with the experimental data of the 
nucleon elastic form factor, one is faced with serious discrepancies, which may 
be due to two important issues, namely the non-relativistic character of the 
model and the smallness of the proton radius. 

In the last decade, new and important data on the proton elastic form factors 
have been obtained by polarization experiments at Jefferson Lab. In fact, if one 
performs elastic scattering of longitudinally polarized electrons from unpolarized 
protons, the recoiling proton has both longitudinal {Pi) and transverse (P*) 
polarization with respect to the momentum transfer in the scattering plane 
m- From the ratio Pt/Pi , one obtains directly the ratio of the proton electric 
and magnetic form factors 


P — 


GW) 

GW) 


^J'P 


Pt Ee + E/ 
Pi 2Mp 



(126) 


where /ip is the proton magnetic moment, is the energy of the incident 
electron, E{, is the energy of the scattered electron, 9^ is the electron scattering 
angle and Mp is the proton mass. 

At variance with the dipole fit, this ratio is deviating strongly from unity 
and it is continuously decreasing towards zero, as shown in Fig. |20[ The mea¬ 
surements have been quite recently extended up to = 8.5 GeV^ |158j . but 
they will be discussed later, in comparison with a reanalysis |159j of the last 
three points of ref. |156j . 

The interesting issue is the existence of a dip in the form factor, as it is sug¬ 
gested by the almost linear decrease of the ratio R. Of course, this unexpected 
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Figure 20: (Color on line) The ratio R = Mp (q 2 ) as a function of the mome- 
tum transfer . The data are taken from |152[ I153L115411155L115611157j . 


result has triggered many theoretical calculations based on quark models at¬ 
tempting to explain the strong deviation of R from unity and discussing the 
possibility of a zero in the electric form factor of the proton. 

It is interesting to remind, at least for historical purposes, that the idea 
of a dip in the electric form factor of the proton has been considered 

already at the early stage of the quark model |160j . The problem was that of 
reconciling the Pauli principle with the completely symmetric character of the 
three-quark wave function in the A resonance: in absence of the colour degrees 
of freedom, the only way out was a completely antisymmetric space function, 
which can be obtained introducing a factor {r \2 — ^ 13 ) (^23 ~ ''’ 21 ) (^31 ~ ^ 32 ) 
the wave function. In general, in presence of a completely antisymmetric space 
wave function, it is easy to show [160j that the charge density in the origin p(0) 
is zero and therefore f dqq^G^^{q) = 0, which means that the form factor 
has a dip with a position strongly dependent on the model form factor. This 
idea has not been further considered since a ground state wave function with 
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nodes is barely acceptable and the introduction of colour made it possible to 
have completely symmetric space functions. 

Actually, there is a calculation, prior to the recent R data, which has pre¬ 
dicted a dip in G^{Q'^). It has been shown |161j that the simple Skyrme soliton 
model, with vector meson corrections, leads to a form factor with a dip at 
~ 1 GeV^ and that, after boosting the initial and final proton state to the 
Breit system, the zero is shifted up to 10 GeV^, a value highly compatible with 
the data of Fig. 20 In later versions of the calculation the zero was pushed to 
~ 16 GeV^ [162]. 


Furthermore, in ref. m. a calculation of the electromagnetic proton form 
factors using a relativistic light-cone constituent quark model has been per¬ 
formed. The plot of R derived from their results exhibits a strong decrease with 
that is due to a zero in the electric form factor at 6 GeV^ m- 
As a final remark, we mention the semiphenomenological fit of the nucleon 
form factors performed in ref. |165j . The parametrized formula contains an in¬ 
trinsic contribution, that now can be ascribed to the quark core, and a vector 
meson cloud. For the charge form factor of the proton, the data were avail¬ 
able only up to 3 GeV^, but if one extrapolates the fitted formula to higher 
momentum transfer, one is faced with the remarkable fact of a dip at abont 9 
GeV^. However, in ref. |165j . the intrinsic (quark core) contribution was not 
considered for the isovector Pauli form factor ■ Its inclusion in has been 
performed in a subsequent paper |166j . but in this case the zero, if present, is 
shifted beyond 10 GeV^. 


Coming back to the hCQM, the ratio R, according to Eq. (124), is identi¬ 


cally 1 in the SU{6) limit and it remains practically 1 also in presence of the 
hyperfine interaction. The hCQM described up to know is nonrelativistic and 
it is worthwhile to analyse the effect of introducing relativity, which, as we shall 
see in the next subsection, is quite beneficial for the description of the decrease 
of the ratio R. 


6.2 Introducing relativity 

In order to describe the electromagnetic form factors both in the elastic and 
inelastic (helicity amplitudes) cases, one has to calculate matrix elements of the 
type 

(127) 

where dt/, 4'/;’ are the initial and final baryon states, respectively and is the 
quark current of Eq. (107). The wave functions of the three quark systems 


are determined in the respective rest frames, but the calculation of the matrix 


elements of Eq. (1271 requires that the quark states be boosted to the Breit 


system. As long as the momentum transfer increases, one expects that the effects 
of the relativistic transformations are more important and therefore relativistic 
corrections must be introduced. In the case of the excitation to the baryon 
resonances, the state B has a mass which can be much greater than the nucleon 
one and in this case the recoil is expected to be less relevant. However, for the 
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Figure 21: The charge (a) and magnetic (b) form factors of the proton. The 
dashed curves are the predictions of the nonrelativistic hCQM, the full curves 
are obtained taking into account the corrections given in Eqs. (1291 and (|130 


The data are from a compilation of ref. m- The figure is taken from ref. jSlj 
(Copyright (2000) by the American Physical Society). 


ground states, nucleon and A in the SU{6) limit, the relativistic corrections are 
unavoidable. 

Relativistic corrections to the matrix elements can be introduced start¬ 
ing from the thee-quark wave function in momentum space in the rest frame 
ip{pp,p\) and applying the Lorentz boosts [30] 

3 

^7 = u{p^) 'tjji{pp,px), (128) 

i=l 

where Bi{i = 1,2,3) are the usual Dirac boost operators that transform the 
quark spinors u{pi) from the nucleon rest frame to the Breit one; the quark 
momenta Pi are in the baryon rest frame and the variables Pp,px are conjugate 
to the standard coordinates p, A, always in the rest frame. The hnal state is 
written in a similar way. 

In order to get the first order corrections to the nonrelativistic approach, the 
matrix elements are expanded keeping the first order in the quark momenta, but 
any order in the momentum transfer m- 

This approach can be applied to the calculation of the elastic form factors 
of the nucleon, obtaining simple analytical forms m 


Ge{Q^) 

/Vr2 

(129) 

Gm(Q") 


(130) 
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Figure 22: The ratio R predicted by the hCQM with the relativistic corrections 
[MZlEI]; the data are from ref. usa. 


where are the electric and magnetic form factors predicted by the 

nonrelativistic hCQM, M is the nucleon mass, E its energy in the Breit frame 
and Je , /m are known kinematic factors [5D] : 






(131) 


where 


^ E 12' 


ti = 


Mm 


9a = oVi + (132) 


6E 


with 


2 I \ll2 |'1QQ^ 

Vi = {m + ) ' ■ (133) 

One sees from Eqs. (|129|) and (|130|) that, besides the multiplicative factors, the 


® = <”* + ■ 


effect of relativity is to rescale the argument of the nonrelativistic form factors. 
It is just thanks to this mechanisms that in ref. USD, the zero of the soliton 
form factor has been shifted to higher values. 

In Fig. (a) and (b) the results for the electric and magnetic form factor 
of the proton are shown. The predictions of the hCQM are higher than the 
data, mainly because the proton radius is about 0.5 fm, however the relativistic 
corrections are beneficial, although there is still a discrepancy with respect to 
the data. 


The interes ting point in connection with Eqs. (1291 and (130) is that the 
ratio R of Eq. (126) acquires a dependence because of the factor [i 


In 
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Figure 23: The transverse helicity amplitudes for the D13(1520) resonance. 
The full curves are the predictions of the hCQM, reported in [45l|46]. The dot- 
dashed (dashed) curves are the results with the relativistic corrections in the 
Breit (equal velocity) frame (From ref. |147j l. 


Fig. 1^ the ratio R obtained applying the boosts to the nonrelativistic hCQM 
is given by the full curve, which deviates from the horizontal line, that is the 
prediction of the nonrelativistic hCQM. 

One can conclude from Fig. that the decrease of R can be ascribed to 
relativistic effects mim- This can be understood if one considers that the 
nonrelativistic results for the form factors are obtained in the nucleon rest frame 
and that the charge and magnetization densities, which are equal in the rest 
frame, behave in a different manner when a Lorentz boost to the Breit frame is 
applied. 

A similar approach can be applied also to the helicity amplitudes predicted 
by the hCQM. In this case one obtains again simple analytical forms |147] . 
which are omitted here for simplicity. The corrections have been calculated for 
the transverse excitation of various resonances, both in the Breit and in the equal 
velocity frames, showing that in any case the deviation from the nonrelativistic 
predictions are very small. As an example, we show in Fig. the results for 
the transverse helicity amplitudes of the D13(1520) resonance |147j . A similar 
situation is valid also for the 811(1535), 811(1650), 831(1620) and D33(1700) 
states 
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7 The relativistic hCQM and the elastic form 
factors 

7.1 The relativistic formulation 

The inclusion of relativistic effects applying Lorentz boosts to the nonrelativistic 
wave functions is an approximation which is beneficial, as shown in the preceding 
Section, but certainly not sufficient. What is needed is a relativistically covariant 
theory, which can be achieved by means of the Bethe-Salpeter approach or 
considering one of the forms of relativistic dynamics introduced by Dirac |lb8] . 

A relativistic covariant theory is characterized by its behaviour under the 
transformations of the Poincare group. To this end it is sufficient to consider the 
generators of the infinitesimal transformations = 0,1,2,3) for the space- 
time translations and Ji and Ki{i = 1,2,3) for the space rotations and boosts, 
respectively. These generators obey to the following commutation relations 
(Po = H) 


[J„P,] = [J„il]=0, (134) 

[Tj, Tj] = (1^^) 

[K,,Pj]=5,,H, [K,,H]=iP,, [K,,Jk]=e^kjK^. (136) 

The problem of building a relativistic theory is equivalent to Hnding a solu¬ 
tion to the Eqs. (134), (135) and (136). Some of the ten quantities P^, Ji,Ki are 
complicated by their dependence on the hamiltonian, while the remaining ones 
are interaction free. As shown in ref. |168j . there are three different types of 
solutions, called instant, point and front form, differing in the type and numbers 
of interaction free quantities. In practical calculations all of them can be used, 
but here we shall concentrate ourselves on the Point Form (PF), for which the 
angular momentum operators Ji and the boosts Ki are interaction free |169] . 
while the four-momentum operator contain the interaction. 

The general three quark states are defined on the product space 
of the one-particle spin-1/2, positive energy representations Hi = L^{R^) x 
(i=l,2,3) of the Poincare group |1701 [3^ and can be written as 

|Pl)P2,P3) Al, A 2 , A 3 ) = |pi, Ai)|p 2 , A2)|P3) A 3 ), (137) 


where pi is the four-momentum of the i-th quark and Xi the z-projection of 
its spin. The states \pi,Xi), i=l,2,3 are given by the Dirac spinors u(pi). We 
introduce the so called velocity states ITHmiTTO) as 

|u, fcl, fc2, fcsj Ai, A2, A3) = Pb(,j) |/ci, ^2, ^3; Ai, A2, A3 )o, ( 138 ) 

where the suffix 0 means that the quark three-momenta ki satisfy the rest frame 
condition ki = 0. With canonical boosts Pb(„), the transformed quark 

tetramomenta are given by pi = UB(v)ki and satisfy the relation 




pfl 

~M 


H e(fci) 


2=1 


(139) 
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where e{ki) is the rest frame quark energy, is the observed nucleon 4- 
momentum and M its mass. In this way the conditions for the PF are satisfied 
[T691I53]. so that the rest frame quark momenta ki and the quark spins undergo 
the same Wigner rotation when the boost is applied. A particularly appealing 
consequence is that the combination of the angular momentum states can be 
performed using the Clebsch-Gordan coefficients as in the nonrelativistic theory 

ITMlfTTTI . 


The hCQM can now be formulated in a consistent relativistic approach. The 
model hamiltonian of Eq. (82) is substituted with the mass operator M [32] 


M 


\Jkf + m'^ 

i=l 


T 

X 


“t“ GX “t“ Hhyp^ 


(140) 


the rest frame quark momenta ki satisfy the condition X]i=i ~ In. the 
calculation, the internal quark momenta are substituted with the Jacobi ones 
Pp = :^iki — k- 2 ), p\ = + k 2 — 2 /c 3 ), which are compatible with the 

rest frame condition and undergo the same Wigner rotations as the internal 
momenta ki. 

The mass operator M can be considered as the result of a Bakamjian-Thomas 
(BT) construction |172j . in which the interaction is introduced by adding to the 

free mass operator Mg = + ttP an interaction Mj in such a way 

that Ml commutes with the non interacting Lorentz generators Jk,Ki. The 
momentum operator Pp can be written as MVp, where the velocity operator 
Vp commutes with the total mass operator M, which is independent of the non 
interacting four velocity |173j . 

The mass operator M can be diagonalized and the resulting baryon spectrum 
turns out to be not much different from the nonrelativistic one, with only a 
slight modihcation of the fitted parameters a and r. One obtains in this way 
the three-quark wave functions tp{kp,k\). The bound states in the rest frame 
are 

tlj{kp, kx)\ki, k2,k3, Ai, A 2 , A3)o, (141) 


which are then boosted to any reference frame according to the Eq. (138) 


UB(v)'4’ikp, kx)\ki,k 2 , fcs, Ai, A2, A3)o = ipikp, kx)\v, fci, fe, ^3, Ai, A2, A3), 

(142) 

the baryon states in any frame are then given by superpositions of velocity 
states. 


7.2 The elastic form factors in the relativistic theory 

In order to describe the electromagnetic properties of the baryons, it is necessary 
to introduce the quark current. For the i-th quark the current is given by 

u(,Pi)ji,{i)u{p'i) = u{pi)e^p{i)u{p'i), (143) 
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{u{pi) are the quark spinors), with which one can calculate the matrix elements 
of the total quark current in the space of the single quark free spinor states: 


{Pi,P2,P3, Ai, A2, A3I J^|pi,P2,P3, a;, A2, A3) 
'EiU{p^)jtl{^)uip'i)u{pj)u{p'^)S{pJ - p'^)u{pk)u{p'^^)S{pk - p'k), 


(144) 


the current conservation is ensured by applying the simple substitution J'^ = 
Jfj. — where is the photon virtual momentum. 

The elastic form factors are extracted from the matrix elements of the current 
between two nucleon states of the type of Eq. (1421, provided that the boosts 
are chosen in such a way to transform the initial and final nucleon states to the 
Breit system. In this way one can calculate the predicted values of the nucleon 
form factors in the hCQM in the relativistic version [35] . The theoretical results 
are compared with the experimental data in Fig. |24| 

The ratio R turns out to deviate from unity more strongly than in the 
semirelativistic analysis shown in Sec. 6.2, in which it becomes flat at a value 
of about 0.6. 

The relativistic calculations are not very far from the experimental data and 
this is a nice achievement if one thinks that the curves in Fig. 24 are predictions, 
however the residual discrepancy indicates that something is missing in the 
theoretical description. 

We should remind that the CQM accounts only for the contribution of a 
quark core and that the proton radius corresponding to the hCQM wave func¬ 
tions is much lower that the experimental one. In fact, according to the analysis 
of the helicity amplitudes for the electroexcitation of the nucleon resonances, 
the comparison with the experimental data and the evaluation of the pion cloud 
terms has shown that there is the need of a meson or quark-antiquark cloud 
surrounding the three valence quarks. In this respect it is interesting to note 
that the fit performed in ref. |I65j is based on expressions of the form factors 
including both an intrinsic and meson contribution; moreover the intrinsic con¬ 
tribution corresponds to a radius ranging from 0.34 fm to 0.55, according to 
the type of form factor used, dipole or monopole, respectively. There are var¬ 
ious attempts in the literature to add meson contributions to the quark core 
ones, but the problem of a consistent treatment of the quark-antiquark pair 
creation mechanism, that is of unquenching the CQM, is very complicated and 
only recently there has been an important improvement for the case of baryons 

[TTmrrnsiri^ . 


There is another relevant feature which is missing. The valence quarks can 
be viewed as effective degrees of freedom which take into account implicitly 
complicated interactions involving also gluons and quark-antiquark pairs. As 
a consequence of such interactions, the quarks can acquire a mass and also a 
size. This statement is supported by the analysis of the deep inelastic electron- 
proton scattering performed in ref. [I78j within the Bloom-Gilman duality: the 
inelastic scattering is due an elastic scattering on constituent quarks with an 
approximate scaling rule given by the quark form factor. Therefore it seems 
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Figure 24: The nucleon elastic form factors predicted by the relativistic hCQM 
(full line) [32]. The data are taken from the reanalysis made in m of the 
data from m- The G^ data are obtained from the G^ ones and the fit |174j 
of the Jlab data on the ratio R; the data for G^ and G)(f are from refs. |176j 
and m, respectively. The figure is taken from ref. [32| (Copyright (2007) by 
the American Physical Society). 


reasonable to take into account the possibility of quark form factors also in the 
description of the elastic electron-nucleon scattering. 

A detailed expression of these quark form factors is not known, then one 
can assume some phenomenological form and use it in order to get a better fit 
of the elastic nucleon form factors. Presently the role of these form factors is 
actually to parametrize the intrinsic structure of the constituent quarks, but 
also any other effect which is not included in the theory, as for instance the 
quark-antiquark pair or meson cloud contributions. 

Quark form factors as superpositions of monopole and dipole functions has 
been used within the relativistic hCQM [35], leading to a very nice agreement 
with data. The results of the hCQM taking into account the quark form factors 
are shown in Fig. 


25 


where, for the G^,G^ and G^ form factors, the ratio 
to the dipole fit is reported. The fit is performed using the data for the form 
factors G^,G^ and G^ and for the ratio R = /r^G^/G^, while the curve 
for G^ is derived from the ratio R and G^. The resulting ratio R is given in 
Fig.[# 
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Figure 25: The nucleon elastic form factors obtained by the relativistic hCQM 
(full line) with quark form factors [32]. The data are as in Fig. 25 The figure 
is taken from ref. 


(Copyright (2007) by the American Physical Society). 


The theoretical description of the experimental data is quite accurate, spe¬ 
cially considering the fact that any discrepancy is certainly enhanced by the fact 
that the ratio to the dipole form factor is presented. 

7.3 The new data and the problem of a dip in the electric 
form factor of the proton 

Quite recently, the measurement of the ratio R has ben extended up to = 8.5 
GeV^ |158j . The curves of ref. [32] have been extrapolated, without any new 
parameter fit, to = 12 GeV^ and compared with the new data [33]. As 
shown in the left side of Fig. [^ the agreement of the extrapolation with the 
new data is quite good, although some doubt is cast because of the large error in 
the last point. The latter is compatible with a dip, while the theoretical curve 
actually continues to decrease, but there seems to be no indication of a zero at 
finite values of Q^. 

The high behaviour of the proton charge form factor can be studied in 
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Figure 26: The ratio fipG^{Q‘^)/G^{Q^ 

(full line) with quark form factors (32] ■ The data are as in Fig. 20 
is taken from ref. 


obtained by the relativistic hCQM 

The figure 


(Copyright (2007) by the American Physical Society). 


an alternative way by introducing the ratio 


Fp = Q 


23m 


(145) 


where Ff{Q'^) and F 2 {Q'^) are, respectively, the Pauli and Dirac proton form 
factors. According to the analysis of Ref. [ns [Ho], such ratio should reach 
an asymptotic constant value. It is interesting to note that in correspondence 
of a zero of the proton form factor G^((3^), Fp must cross the value 4Mp [55] . 
where Mp is the proton mass. Using the curves for the ratio R and for the 
proton magnetic form factor G^((5^) evaluated in |33|, it is easy to determine 


27 


FfiQ^) and Ff (Q^). The results are shown in the right side of Fig. 

Of course, the two plots in Fig. [27] contain the same information and no 
further comment is needed. 


The new data by ref. m have some problem of compatibility with the 
previous ones. In fact, the trend of the ratio was, mainly thanks to the last 
points of ref. |156) . a linear decrease to a zero given by the fit HHII 


R = 1.0587- 0.14265Q2 


(146) 


As for the data in ref. |158j . the last point, with a very large error, is compatible 
with such trend, but the other two definitely no. This question has been studied 
in ref. m. where the last points of ref. |156j have been reanalyzed, with the 
result that the measured values of R for 3.5 GeV^ < < 5.6 GeV^ should be 
increased, improving the consistency with the higher data. The presence of 
a dip in the charge form factor of the proton is still an open problem, which 
will hopefully solved by the planned experiments at the 12 GeV upgrade of the 
Jefferson Lab electron accelerator. 
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Figure 27: Left: The ratio fipG^{Q‘^)/G^{Q‘^) obtained by the relativistic 
hCQM (full line) with quark form factors [33]. The data are as in Fig. 20 and 
from |158| . Right: the ratio Fp = (Q^)/Ff (Q^) obtained by the relativistic 

hCQM (full line) with quark form factors [33]; the data are as in the left side. 
The figure is taken from ref. [33] (Copyright (2010) by the American Physical 
Society) 


Coming back to the hCQM, the theoretical curve reported in Fig. 27 shows 
no evidence of dip, in agreement with the new data and the modified [156] 
points. It should be noted that in ref. |158j the present data have been fitted 
by a curve for which the dip, if any, is pushed well beyond 10 GeV^. 

The problem of the R ratio has been dealt with also by other groups involved 
in the building of CQMs. 

The Rome group [28] have calculated the R ratio using the model of ref. [69] 
in the front from dynamics, obtaining a decrease with given by the Melosh 
rotations. Using the GBE model with point form dynamics, the Pavia-Graz 
group [M] [3S] reproduced the form factor up to 4 GeV^ and a decrease of the 
ratio. A similar behaviour |36| is obtained by the Bonn group within the Bethe- 
Salpeter approach with instanton interaction. For further details, the reader is 
referred to refs. [18211183] . 


7.4 Relativistic transition form factors 

As quoted in Sec. 6.2, the helicity amplitudes seem to be practically unaffected 
by the application of Lorentz boosts, with the possible ecception of the N — A 
excitation. These statements can now be checked by calculating the helicity 
amplitudes with the relativistic hCQM described in the previous sections. These 
work is presently in progress and for the moment there are only preliminary 
results on the N — A transition, which is expected to be more sensitive to 
relativistic corrections, as it happens with the elastic form factor. One should 
not forget that in the SU{6) limit, the nucleon and the A states are degenerate 
and belong to the lowest shell. 

The procedure for the calculation of the N — A helicity amplitude is the 
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Figure 28: The transverse helicity amplitudes for the TV — A excitation. The 
results of the relativistic hCQM (full curves) are compared with the nonrelativis- 
tic ones (dashed curves) and with the experimental data, a) ^f /2 amplitude; b) 
amplitude. The data are from [1311 ll21j ('From ref. [47]). 


same as in the case of the elastic form factor, apart from the fact that in this 
case the final baryon state is given by a A. The results |47| are given in Fig. 

As expected, the calculations in the fully consistent relativistic hCQM pro¬ 
duce a sensible improvement of the results. There is still a lack of strength at 
low which, as mentioned in Sec. 5.3, is due to the missing quark-antiquark 
pair creation mechanisms. 

For the remaining resonances, the approximate inclusion of relativistic effects 
performed in Sec. 6.2 is certainly preliminary and a program of a systematic rel¬ 
ativistic calculation of both the longitudinal and transverse helicity amplitudes 
is desirable. In fact there is an important experimental program at the 12 GeV 
upgrade of Jlab and the availability of a consistent quark model will be useful 
for the analysis of the high behaviour of inelastic processes, which are ex¬ 
pected to provide valuable information on the short distance properties of the 
nucleon [11511116) . 
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8 Summary and conclusions 


The hCQM provides a good basis for the description of many baryon properties, 
such as the spectrum, the transition amplitudes for the photo- and electro¬ 
excitation of the resonances and the elastic nucleon form factors. 


The model potential (see Eq. (82)) has, as it happens in all CQMs, a domi¬ 
nant iS'{7(6)-invariant part and a perturbative term producing the SU{6) viola¬ 
tion, a structure which is present already in the early Lattice QCD calculations 
|59l [60]. The spin-flavour independent potential with its hC and linear con¬ 
finement terms has strong similarities with the Cornell one m and the static 
qq interaction derived from LQCD in the SU{3) limit [35], while the SU{6) 
violating term is based on the One-Gluon-Exchange mechanism [60] . 
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The model has only three free parameters, which are fixed in order to de¬ 
scribe the experimental non strange baryon spectrum, obtaining a good agree¬ 
ment in particular for the negative parity resonances |45j . Once the parameters 
have been determined, all the baryon states can be consistently constructed and 
used for parameter-free calculations. 

Among these there are the predictions of the hCQM for the helicity ampli¬ 
tudes, which are fairly described. The medium-high behaviour is generally 
well reproduced, while there is a lack of strength al low a feature which 
is in common with all CQMs. The results for the 1/2 amplitudes are usually 
better than for the 3/2 ones; this is particularly true for the 811(1535) trans¬ 
verse excitation, whose behaviour has been predicted |45j three years in advance 
with respect to data mu, USD, Ei], [HU- It should be reminded that the 
calculated proton radius turns out to be 0.48 fm; such a value was obtained by 
the fit of the D13 photocoupling [21], while in the hCQM case is a consequence 
of the fit to the spectrum. 

It should be observed that the hC term 1/x jH [Ml lU plays a crucial role 
in obtaining a good agreement with the experimental data. First of all, there 
is a perfect degeneracy between the first negative parity level 1]^ and the first 
hyperradial excitation Oj, a feature which is only slightly modified by the con¬ 
finement term and is at variance with any two-body potential. This fact provides 
a better starting point for the description of the Roper resonance, although its 
precise position can be given only if flavour dependent interactions are consid¬ 
ered. Moreover, the particular structure of the 1/x spectrum allows to include 
within the first three shells all the observed 4- and 3- star negative parity res¬ 
onances, while in the h.o. case one should involve the 37iw leveles, the more 
so if one considers also the 2-star states urn). Finally, the presence of the hC 
potential leads to more realistic quark wave functions, which allow a better de¬ 
scription of the helicity amplitudes [13 [H]. This statement is supported by 
the fact that the analytical model [93l [94], where the confinement is treated 
as a perturbation, makes use only of the hC wave functions and the resulting 
helicity amplitudes are very similar to those predicted by the complete model 
of Eq. 

The nonrelativistic hCQM has been extended in order to include also isospin 
dependent interactions m and to describe strange baryon resonances m- 

The above results have been obtained in the nonrelativistic version of the 
model and are quite satisfactory, apart from the elastic nucleon form factor. In 
fact, as long as the transfer momentum is increased, the theoretical predic¬ 
tions deviate from the experimental data, both because of the small theoretical 
proton radius and because of an incorrect account of the recoil in passing from 
the rest frame, where the calculations are performed, to the Breit frame which 
is more appropriate for the interpretation of the form factors as related to the 
charge and magnetic distributions. In principle, these considerations are valid 
also for the inelastic transitions to the resonances, but in this case, because 
of the large masses of the resonances, the recoil is less important, apart from 
the lightest resonance, that is the A(1232). On the other hand, one should 
not forget that the nucleon and the A(1232) are, in the SU{6) limit, both in 
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the ground state configuration. Actually, the relativistic effects turn out to be 
important for both the elastic nucleon form factors and for the electromagnetic 
excitation of the A(1232) resonance. 

Relativity can be introduced in the hCQM calculations of the form factors, 
both elastic and inelastic, by simply applying Lorentz boost to the nonrelativis- 
tic three quark states [5DI fSTl ITT7] . As expected, the helicity amplitudes for the 
higher resonances are not much affected by this relativistic corrections. On the 
contrary, the boosted elastic nucleon form factors are improved but are still not 
in agreement with the experimental data. 

A consistent way of constructing a relativistic quark model is provided by 
the forms of relativistic dynamics introduced by Dirac |168j . We have chosen 
the Point Form (PF), which allows to combine the angular momentum states 
in the same way as in the nonrelativistic case and considered a mass operator 
according to a Bakamjian-Thomas (BT) construction |172j and containing the 
quark interaction of the hCQM, that is a hC potential and a linear confinement. 
Refitting the parameters, which turn out to be very similar to the nonrelativistic 
ones, it is possible to predict the elastic form factors [5^ I33j . which are now 
very near to the data. Also the predictions for the A excitation are improved if 
the relativistic PF is used. A very good description of the elastic form factors 
is achieved if quark form factors are introduced |32| . 

An important issue in connection with the elastic nucleon form factors is 

qP ('Q2'\ 

the observed behaviour of the ratio R = iip (q 2 ) with increasing values of 
the momentum transfer. In the nonrelativistic model, this ratio is practically 
one, while already the simple application of Lorentz boosts leads to a ratio 
R which deviates from one, agreement with the idea that the depletion is a 
relativistic effect |167L I31j . In fact the relativistic hCQM, with the inclusion 
of quark form factors, is able to describe the behaviour of R as well [32]. As 
for the possible presence of a dip in the electric form factor, the situation is 
not yet clear and hopefully it will be settled by the future Jlab experiments. 
Extending the relativistic hCQM calculations up to 12 GeV^, there seems to be 
no indication of a zero in the proton electric form factor. 

The hCQM shares with all other models the drawback that the predicted 
resonance states have zero width, that is there is no coupling to the continuum. 
To this end one can introduce some mechanism for the production of mesons 
and in this way it is possible to calculate the baryon decays as for instance in 
refs. |22l|50l|5ll|52l|53|. However the coupling to the continuum is important 
also for its consequences on the spectrum, in particular for the possible pres¬ 
ence of meson cloud or quark-antiquark components in the states. The latter 
feature can be introduced by simply considering higher Fock components in the 
baryon states, either in the form of a meson-baryon or Aqq configurations. A 
more consistent procedure is preferable and this can be achieved by general¬ 
izing the quark interaction in order to take into account the coupling to the 
continuum and then perform ab initio calculations of the various quantities of 
interest. An attempt in this direction is the work of ref. |106j . where however 
quark and meson-baryon degrees of freedom are treated on the same footing. 
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An important progress in this direction is provided by the recent unquenching 
of the CQM formulated in refs. |1()7L 11081 llOhj . in which in particular the quark- 
antiquark pair contributions to the spin and the flavor asymmetry in the proton 
have been calculated consistently. As already mentioned, the quark-antiquark 
pair creation mechanism is expected to be the main contribution to the elec¬ 
tromagnetic excitation strength at low Q^. In any case, these works open the 
way to a consistent quark model calculation of the baryon spectrum, the strong 
decays, the meson electroproduction and the elastic form factors. 
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Appendix A. The hyperradial wave functions 

The hyperradial wave functions (x) are solutions of the hyperradial Schodinger 

equation (74) of the text 


dx"^ ' X dx 


- — - ^ ^^ ] = - 2m [E - V 3 q{x)] iIjji,{x). (1) 


they are normalized according to the condition 

nOO 

/ dX X^\'tjj.y„{x)\‘^ = 1 . 

Jo 


( 2 ) 
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Table 1: The hyperradial wave functions 'ip-y^{x) for the harmonic (h.o.) and 
hyperCoulomb (hC) potentials, up to the second energy shell. In the h.o. case 
a = ^/SKm/h, Nh.o. = 2:/ + 7 and an overall common factor is omitted. For 
the hC case Nhc = 1 ^ + 7 and the constant g is given by rm/{v + j + 5/2). 



The reduced wave function is defined as 

-tp^„{x) = 

and obeys to the reduced equation 


Uj^(x) 

x5/2 


( 3 ) 


cPu^^(x) , , 15,u^^(x) 2m / m / ^ /-n 

-^^-bh + ^)-^]^^ + -^[E-V3gix)] u^^x) = 0; (4) 

the behaviour for small values of x is given by 

u^v{x) ~ . (5) 

In the case of the h.o. interaction V 3 q{x) = Zjlkx^, the asymptotic behaviour 
is determined by the gaussian factor 


U~^u{x) 


( 6 ) 


for the reduced wave function one can then assume the form, apart from a 
normalization factor. 


-^ 71 / 


{x) = 


— T,5/2+7„-a^a:' 




( 7 ) 
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( 8 ) 


In order that Eq. Q be satisfied, the function P must obey to 

,a ^.dP P 

+ (- 1)-^ - \-V— — 0, 

y dy y 


where the variable y = has been introduced and 

_ 3 + 7 e 
^ 4 ^’ 


a = 3 + 7 , 


( 9 ) 


with € = ^E. Eq. (j^ is the hypergeometric confluent equation and P is then 
given by the hypergeometric confluent function + 7 ,y). In order to 

preserve the correct behaviour for small and large x, ly has to be a non negative 
integer and P(y) is a polynomial of order v. The energy values are then given 

by 

E = (3 + 2i/+ 7 )?iw, (10) 

with w = 

The h.o. wave functions for the first three shells are reported in Table [T] 
For the hC potential — ^, one can proceed in an analogous way. The reduced 
equation @ is written as 

(x) 15 (a;) 2m 

-[ 7(7 + 4 )-—] ^ [E-V-iq{x)\ u^^{x) = 0 . ( 11 ) 

In this case the asymptotic behaviour is determined by the factor where 

= — V2mE, ( 12 ) 

(h is taken equal to 1). One can then assume 

u-yiy{x) = e~^^P{x). (13) 

Eq. Q is satisfied if the function P obeys to 


d^P ^,dP P 

^ + (--l)^+z/- = 0, 


dy 


y 


dy 


y 


(14) 


where the variable y is now defined as y= 2 gx and the parameters a and v are 
given by 

r r, 5 mr 

a = 5 + 27 , —v= -+7 -. 

2 g 

From the last formula one obtains finally [94] 

mr 

9 = 5- 


+7 + 


(15) 

(16) 


and 


E = - 


2m 


2[f+ (7 + ^)(7 + ^ + 5)]’ 


(17) 


The quantities P(y) are associated Laguerre polynomials (see Eq. (801 in the 
text). Also for the hC wave functions for the first three shells are reported in 
Tabled! 
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Appendix B. The baryon states 

The baryon states are superpositions of 5'C/(6)—configurations, which can be 
factorized as follows (see Eq. Q of the text) : 

^colour ‘ \spin ' ^isospin ^ 35 (P,A). ( 1 ) 

As already mentioned in the text, the various parts must be combined in 
order to have a completely antisymmetric three-quark wave function. To this 
end it is necessary to study the behaviour of the different factors with respect 
to the permutations of three objects (that is with respect to the group S 3 ). In 
general, any three particle wave function belongs to one of the following symme¬ 
try types: antisymmetry (A), symmetry (S), mixed symmetry with symmetric 
pair (MS) and mixed symmetry with antisymmetric pair (MA). 

For the colour part Ocoiour one must choose the antisymmetric colour singlet 
combination. 

The three-quark spin states are dehned as: 


Xms 

= I((2.2)1>2^2^’ 

(2) 


1, 1,1, 


Xma 


(3) 


1, 1,3, 


Xs 

= I((2>2)^’2^2^- 

(4) 


The antisymmetric combination is absent because there are only two states 
at disposal for three particles. 

Similarly one can define the isospin states ipMS, (I’mA, 4>S- 
If the interaction is spin and isospin (flavour) independent, one has to intro¬ 
duce products of x~ and </>— states with definite S 3 — symmetry. Here we give 
the explicit forms only for the case that both factors have mixed symmetry, the 
remaining ones being trivial: 


[XMA(t>MA + XMS(t>Ms], 

^MS = [xma4>ma - Xms4>ms], 

^MA = \xma4>ms + Xms^’ma], 

^A = ^ [xma4>ms — Xms^'ma]- 
The space wave function is given by 

'4^3q{PT^) — (Up, Haj C)) 


(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 


77 


where ^ = 2n + Ip + ly, the hyperspherical functions are given by [55] 


— '^Ipiripi^p) yixm^i^x) ^ ^-P/v ^(?)) 


( 10 ) 


where 


with 


= Cr,i^ix{cosO^^{sinS,p 


Cnlplx — 


/2(27 + 2)r(7 + 2 -n)r(n + 1 ) 
r(n + ;p + 3/2)r(n + ZA + 3/2 


( 12 ) 


and P^'^{z) is a Jacobi polynomial. 

The symmetry property of the wave function Eq. is determined by the 
hyperspherical part only, since the hyperradius x is completely symmetric. In 
Table we report the combinations of the hyperspherical harmonics having 
definite S'3—symmetry. 

In Tables and we give the explicit form of the three-quark states with 
positive and negative parity, respectively. In these Tables the hyperradial wave 
functions are solutions of the hyperradial equation Eq. (74) of the text; 
their form depends of course on the hypercentral potential. 
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Table 1: Combinations (y[ 7 ]Zp;;,)s 3 [M] of the hyperspherical harmonics 
that have definite 5 ' 3 —symmetry. For simplicity of notation, in the third column 
we have omitted the coupling of Ip and lx to the total orbital angular momentum 
L. Each combination is labelled as specifying the total orbital angular 
momentum L, the parity P and the S' 3 —symmetry type t = A, MA^ MS, S. 


7 

Lf 

(^[7]Cc)s3 

^3 

0 


^[0]00 

S 

1 


y[l]10 

MA 



y[i]oi 

MS 

2 

2+ 

^[^[2]20 + ^[2]02] 

S 



^[2] 11 

MA 



^[^[2]20 - ^[2]02] 

MA 



y[2]ll 

A 


0^/ 

^[2] 11 

MA 



^[2]00 

MA 
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Table 2: Three-quark states with positive parity [94]. The second, third and 
fourth columns show the angular momentum, parity and S' 3 -symmetry, , the 
spin, S, and isospin, T. States are shown in the last column and are written in 
terms of the hyperradial wave functions, of hyperspherical harmonics, 
(Y[..y])s 3 of Table and of the spin and isospin states. 


State 


s 

T 


SU( 6 ) configurations 

Pll 

Os 

1 

2 

1 

2 


V'oo y[o]oo f^S 


oj 

1 

2 

1 

2 


V'lo b^[o]oo f^s 


oj 

1 

2 

1 

2 


■020 ^[0)00 f^S 



1 

2 

1 

2 


V'22 [^[2]00 f^MS + b]2]ll ^^ma] 



3 

2 

1 

2 

V'22 

(^[2]2o ~ ^^2102) 0 MS + ^[2]ii 0 ma] Xs 

P13 

2+ 

1 

2 

1 

2 

"022 (^[2]20 ~ ^[2]02) f^MS +b]2]ll ^ Ma \ 



3 

2 

1 

2 

V'22 

77f [75 (^[2]20 - ^[2]02) 0 MS + b"[2]ii 0 Ma] XS 


Om 

3 

2 

1 

2 


022 ^ [^[ 2]00 0 MS + ^[2]11 0 Ma] XS 


2+ 

1 

2 

1 

2 


022 [b^[2]20 + ^[2]02] f^S 

P15 

2+ 

1 

2 

1 

2 

7^ [772 (^[2]20 ~ ^^2)02) f^MS + ^[2]11 ^ Ma \ 


2+ 

3 

2 

1 

2 

V'22 

^ [:^ (^^2)20 ~ ^^2102) 0 MS + ^[2]ii 0 ma] Xs 


2+ 

1 

2 

1 

2 


022 [^[2]20 + ^[2]02] 

P17 

2+ 

3 

2 

1 

2 

V'22 

77 f [:^ (^[ 2]20 - ^[2]02) 0 MS + ^[2)11 0 ma] XS 

P31 

2+ 

3 

2 

3 

2 


022 ^ [(^[2120 + ^[2]02] XS 0 S 



1 

2 

3 

2 


022 ^ [^[2]00 XMS + b"[2]ll XMa] 0 S 

P33 

Oj 

3 

2 

3 

2 


000 b]o]oo xs 0s 


oj 

3 

2 

3 

2 


010 ^[0)00 xs 0s 


oj 

3 

2 

3 

2 


020 ^[ 0)00 Xs 0s 


2 s 

3 

2 

3 

2 


022 ^ [^[2]20 + ^[2]02] XS 0 S 


2 + 

1 

2 

3 

2 

V'22 

^ (^^2120 ~ ^[ 2 ) 02 ) XMS + ^[2]11 XMa] 0 S 

P35 

2 + 

1 

2 

3 

2 

V'22 

(^[2]20 ~ ^[2)02) XMS + ^[2]11 Xma] 0 S 


2+ 

3 

2 

3 

2 


022 [b^[2]20 + ^[2]02] XS 0 S 

P37 

2 + 

3 

2 

3 

2 


022 [^[2]20 + ^[2]02] XS 0 S 
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Table 3 : Three quark states with negative parity [ 13 ]. Notation as in Table 


Resonances 


S 

T 




States 


511 


1 

2 

1 

2 

V'l 

1 ^ [^[l]10 + h"[l]01 ^Ms] 



1 

2 

1 

2 

i’2 

1 ^ [^[1]10 + RjlJOl] ^MS 



3 

2 

1 

2 

^11 

1 

V2 

[^[1]10 

i>MA + ^[1]01 

(jiMs] Xs 



3 

2 

1 

2 

tp21 

1 

\/2 

[^[1]10 

ilMA + ^[1]01 

</'Ms] XS 

D 13 


1 

2 

1 

2 

i’ll [^[l]10 + h]l]01 ^Ms] 



1 

2 

1 

2 

V'21 [^[1]10 + ^[l]0l] ^MS 



3 

2 

1 

2 

^/’ll 

1 

V 2 

[>^[1]10 

ilMA + ^[1]01 

(jlMs] xs 



3 

2 

1 

2 

i’21 

1 

V 2 

[^[1]10 

(I’M A + ^[1]01 

(iMs] xs 

D 15 


3 

2 

1 

2 

i’ll 

1 

V 2 

[^[1]10 

4 >MA + ^[1]01 

(^Ms] Xs 



3 

2 

1 

2 

i’21 

1 

V 2 

[^[1]10 

(I’M A + ^[1]01 

</'Ms] xs 

531 


1 

2 

3 

2 

i’ll 

1 

V2 

[>[1]10 

Xma + ^[1]01 

Xms] (is 



1 

2 

3 

2 

i ’21 

1 

V2 

[>[1]10 

Xma + ^[1]01 

Xms] (is 

533 


1 

2 

3 

2 

i’ll 

1 

V 2 

[>[1]10 

Xma + ^[i]oi 

Xms] (is 



1 

2 

3 

2 

i ’21 

1 

J 2 

[>[1]10 

Xma + ^[i]oi 

Xms] (is 
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